TTP12-032 



Analysis of the Higgs potentials for two doublets and a singlet. 

G. Chalons, F. Domingo 

Institut fiir Theoretische Teilchenphysik, 
Karlsruhe Institute of Technology, Universitdt Karlsruhe 
Engesserstrafie 7, 76128 Karlsruhe, Germany 



Abstract 

Wc consider the most general CP-conserving renormalizable effective scalar potential involving 
two doublets plus one singlet Higgs and satisfying the electroweak gauge symmetry. After deriving 
the electroweak-symmetry breaking conditions, we focus on special cases, characterized by specific 
symmetry properties and/or relations to supersymmetry-inspired extensions of the Standard Model 
(e.g. n/NMSSM, UMSSM). We then investigate the question of the reconstruction of the potential 
parameters from the Higgs masses and mixing angles and show that in some specific cases, such 
as the one of an underlying NMSSM, an accuracy at the order of leading-logarithms is achievable 
with minimal effort. We finally study a few phenomenological consequences for this latter model. 
More specifically, we consider how our parameter reconstruction modifies the outcome of two publicly 
available codes : micrOMEGAs and NMSSMTools. We observed noteworthy effects in regions of parameter 
space where Higgs-to-Higgs decays are relevant, impacting the collider searches for light Higgs states 
and the prediction of the Dark-Matter relic density. 



Introduction 



The origin of ElectroWeak Symmetry Breaking (EWSB) stands as one of the critical questions in high- 
energy physics and a central goal of the Large Hadron Collider (LHC) is to reveal its nature. The 
recent discovery of a new massive boson around 125 GeV [l], reported by both the ATLAS and CMS 
collaborations [2j , and supported by the broad excess seen at TeVatron , represents a first step towards 
the identification of the Higgs boson and the measurement of the underlying Higgs potential, a task which 
however only the next generation of colliders will probably complete. Although essentially compatible 
with the Higgs boson of the Standard Model (SM), this new state may already be hinting towards some 
new physics, in that the peaks of the diphoton and ZZ — )■ 4Z decays differ from what one would expect 
in the SM. The stronger signal in the i7 — )• 77 channel, in particular, seems of importance because this 
loop-induced process is particularly sensitive to physics beyond the SM. One should also consider the 
non-observation of events at CMS - although supported by very little statistics - in the iJ — )• rr channel. 
Testing the SM-nature of this would-be Higgs state, inspecting possible deviations in its coupling to SM 
particles shall represent a major undertaking of modern particle physics and a probable probe into the 
mechanism of EWSB. 

The 'Higgs mechanism' involving scalar elementary fields, is the most efficient way to generate masses 
for the fermions and gauge-bosons. Its implementation within the SM is the minimal one: only one 
scalar field, transforming as a doublet under SU{2)l, is introduced to break the electroweak (EW) 
symmetry through its vacuum expectation value (v.e.v.). Nevertheless the Higgs sector is still essentially 
undetermined and there is no reason to stick to minimality if some benefits should emerge from a more 
elaborate scalar sector. For instance, introducing a second Higgs doublet allows for an implementation of 
CP violation through this sector [s): CP violation appears in this context because some of the parameters 
in the potential of the Two Higgs Doublet Model (2HDM) can be chosen complex (non-real). Yet the 
requirements relative to neutral flavor conservation constrain this possibility significantly. Large flavour- 
changing couplings of neutral Higgs bosons can be avoided in the so-called '2HDM of type IV, where 
the Higgs doublets and Hd, of opposite hypercharges Y = ±1, enter separately, and respectively, up- 
and down-type Yukawa terms (at tree level). Another (more exotic) possibility consists in requiring the 
alignement of the Yukawa coupling matrices in flavor space: see fo]. Although such 2HDM's may hold as 
autonomous extensions of the SM, they can also be embedded within more elaborate models: Left-Right 
gauge models and their Grand-Unification Theory (GUT) ramifications - Pati-Salam, 50(10), etc.- offer 
a first framework for this operation, in which the question of CP-violation was originally central j?]. 
From another angle, the well-documented 'Hierarchy Problem' |8| underlines the theoretical difficulties 
for understanding the stability of a Higgs mass at the electroweak scale, with respect to new-physics at 
very-high energies (GUT, Planck scales). Regarding the SM as the low-energy effective theory of some 
more-fundamental model, the quadratic sensitivity of scalar squared masses to new-physics masses would 
lead to a technically unnatural fine-tuning of the Higgs-mass parameter in the more-fundamental theory 
with the radiative corrections resulting from the integrated-out new-physics states. . . Unless new-physics 
appears sufficiently close to the electroweak scale: typically at the TeV scale. Among the proposed 
solutions, Supersymmetry (SUSY) allows to stabilize a scalar Higgs mass at the electroweak scale, due to 
the renormalization properties of supersymmetric theories. However, SUSY being obviously not realized 
in low-energy particle physics, viable SUSY-inspired models need to include SUSY-breaking effects, which 
are parametrized within the Lagrangian through the so-called 'soft terms', generate e.g. mass terms for 
all non-SM particles and trigger the Higgs mechanism. This ad-hoc setup could yet remain an acceptable 
solution to the Hierarchy Problem only if the supersymmetry-breaking scale is near the electroweak 
scale. Other attractive properties of SUSY-inspired models lie in the possibility of one-step unification, 
due to the more-convergent running of SM-gauge couplings in the presence of the enlarged SUSY field- 
content [9], or in the dark-matter (DM) sector, the lightest supersymmetric particle being a stable (or 
long-lived) and viable candidate in the presence of (approximate) R-parity |To|. 

Holomorphicity of the superpotential (cancellation of gauge-anomalies) dictates the requirement for at 
least two SU{2)l Higgs doublets in a SUSY-inspired model, intervening in a Type II 2HDM fashion, so 
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that both up-type and down-type masses be generated. The simplest implementation of a SUSY-inspired 
SM, known as the Minimal Supersymmetric Standard Model (MSSM) confines to this minimal 2HDM 
requirement. There, the quartic Higgs couplings are determined by the EW gauge couplings, which results 
in tight constraints on the tree-level mass of the lightest Higgs boson: the latter is indeed bounded from 
above by the Z'^-boson mass Mz- Radiative corrections improve this feature and can arrange for fairly 
heavy Higgs masses provided the SUSY-scale is large enough, see for example p2] . Yet this last necessity 
tends to conflict with the naturalness-dictated < 1 TeV SUSY-breaking scale. Accommodating for a 



Higgs state at 125 GeV in the MSSM hence severely constrains the parameter space of this model 13 
Another criticism to this minimal setup, the so-called '/i-problem' |14| , points out the necessity of tuning 
a supersymmetric mass-term, the conventionally-baptized /i parameter, at the electroweak/TeV scale in 
order to ensure EWSB: being of supersymmetric origins, this parameter is in principle unrelated to the 
SUSY-breaking scale and would thus coincide with it out of sheer coincidence. 

The introduction of an additional gauge-singlet superfield S addresses both shortcomings of the MSSM. 
The /i-term can indeed be generated effectively through a XSH^ ■ H^i term when the singlet takes a v.e.v. 



s: ^eff. = As 15 . Concerning the lightest Higgs mass, the presence of a new superfield coupling to the 
Higgs doublets induces additional contributions to the Higgs mass matrix, so that the MSSM limit can 
be exceeded, already at tree-level |l6|[T7] . It is also worth to mention that the lightest CP-even Higgs 
state in this context might well be dominantly of a singlet nature, hence, the singlet decoupling from 
SM-fermions and gauge bosons, essentially invisible at colliders: the SM-like Higgs state would then be 
the second lightest and a small mixing effect with the singlet would thus shift its mass towards slightly 
higher values. In short, radiative corrections are no longer the only mechanism able to generate a SM-like 
Higgs-state heavier than Mz in such a singlet-extension. 

The simplest version of such a model with singlet-enlarged superfield content is known as the Next-to- 



Minimal Supersymmetric Standard Model (NMSSM) |18[ 19 . It relies on a Z3 discrete symmetry in 
order to forbid all dimensional parameters (including ji) in the superpotential, so that the soft-terms 
provide the only relevant scale in the scalar potential, triggering the EWSB. Several other SUSY-models 
engaging a singlet in addition to the two Higgs doublets are to be found in the literature, including the 
nearly Minimal Supersymmetric Standard Model (nMSSM, sometimes MNSSM) |20[|2l|, J7(l)'-extended 



MSSM's, with their simplest version known as the UMSSM p2], models based on the Eq exceptional 



group 123], SUSY/compositeness hybrids, such as 'fat Higgs models' 24 or models using the Seiberg 
Duality |25|, etc. 

In the present paper, we aim at studying the effective Higgs potential involving 2-doublet-|-l-singlet Higgs 
fields. The relations between physical input, represented by the mass matrices and mixing angles, and 
the parameters of the potential, as well as with the trilinear Higgs couplings, shall be at the center of 
this discussion, in view of a possible reconstruction of the potential from such input, at, and beyond, 
leading order (LO). Similar analyses for the 1-doublet setup |26j, or the 2-doublet setup, for instance 
in [2^, with the MSSM as a background-model, have already been proposed in the literature. Given 
that the singlet-extensions of the MSSM offer a natural origin to our 2-doublet-|-l-singlet setup, we 
shall refer and return explicitly to such models in the course of our discussion: specific attention will 
be dedicated in particular to the n/NMSSM or the UMSSM. Most of our discussion should however be 
generalizable to other models resulting in a 2-doublet -|-l-singlet Higgs potentiaQ as long as matching 
conditions or /and symmetry properties are satisfied. The first part of the present paper shall be dedicated 
to the presentation of the general framework, including notations, the discussion of residual symmetries 
and the pattern of EWSB leading to the Higgs spectrum. In the second part, we shall focus on the 
question of the reconstruction of the potential from a measurement of the Higgs masses and mixing 
angles: beyond the general case where a large number of undetermined parameters remain, the possibility 
of a reconstruction in constrained models will be discussed at leading order. The analysis of the large 

^We have already referred to Left-Right models and their GUT extensions as an alternative approach to the 2HDM 
framework. Note that the addition of a SM-gauge singlet is essentially an undemanding requirement and may be arranged 
within such models as well. 
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logarithms appearing in the Coleman- Weinberg 28 29] approach shall convince us, in particular, that 



a full reconstruction at the order of leading logarithms should be achievable in the Za-symmetric case 
represented by an underlying NMSSM. Concentrating on the NMSSM in the last part, we shall analyse 
the phenomenological consequences for this model, both in terms of constraints from Higgs-to-Higgs 
decays and computation of the Dark-Matter relic density. The decay 77 [30j will also be revisited, 

although little impact is expected there. This phenomenological analysis will rely on the numerical output 



of several public codes, including NMSSMTools_3 . 2 . HbT 32 , inicrOMEGAs_2.4. 1 (33,34 and a version of 



Sloops 35 36 adapted to the NMSSM 37 



1 Two Higgs doublet plus one singlet potential 
1.1 General parametrization 

New-Physics (NP) effects are most conveniently encoded in terms of effective Lagrangians. Under the 
guidelines of Lorentz and gauge invariance, as well as possible additional symmetries, one can write a list 
of all the operators, classified according to their mass-dimension. For the two SU{2)l doublets and the 
singlet, we shall use the notations (with v^^ Vu and s representing the v.e.v.'s of these fields): 



r;rf + (/i0 + m0)/^/2' 



H,, 



Hi 

+ {hi + 



, 5 = s + (/i° + m°)/\/2 



The most general Higgs potential involving these fields and compatible with the electroweak gauge sym- 
metry then reads, when one restricts to renormalizable terms: 



eflf. 



("^12 



+ ^aIHu ■ Hd\ + 



As 



Hu ■ Hdf' + (Ael-ffnP + X7\Hd\^)Hu ■ Hd + h.c. 



4 



+ 



2 3 3 ' ' 4 

1 2 u TT I \ M o* "2 ZJ U I \ A/ r<2 i 



AudSHu ■ Hd + AudS*Hu • Hd + \m\S\ Hu ■ Hd + Xp S* Hu ■ Hd + Xp S Hu ■ Hd + h.c. 



+ A^l^^i^^l^ + XUSl'^m'' + {AusS + X}>S')\Hu\' + {AdsS + X'i,S^)\Hd\' + h.a 



[1.2) 



The first two lines comprise the usual 2HDM potential, the third one, the pure-singlet terms and the 
latter two, the singlet-doublet mixing-terms, rn^^, fn^j^, Ai, A2, A3, A4, m|, k? , Ap and Ap are (10) real 

parameters, while mfa, A5, Ae, A7, At, A5, A5, k|, k|, Aud, Aud, Am, A^^, A^^, Aus, Ads, Ap and A^ 
are (19) in-principle-complex parameters. One parameter (e.g. At) is superfluous and may be absorbed 
in a translation of the singlet; three others (m|, rn'jj^ and rriff^) can be traded for the field vacuum 
expectation values through the minimization conditions. From now on, we will consider, for simplicity, 
that all the parameters are real, hence barring the possibility of CP-violation. (We will however continue 
to refer to the 19 potentially non-real parameters as 'complex' parameters.) 



1.2 Symmetry classification 



By imposing additional symmetries, the form of the potential in Eq.(1.2) can be further constrained 
at the classical level and the remaining parameter^ A^^' will be called 'classical' parameters. At the 
quantum level, all the eliminated terms A^™ may reappear, in principle, if the symmetry is broken, either 
directly by the quantum fluctuations, or spontaneously, when the fields acquire v.e.v.'s. In the later 



We shall use the notation 'Ai' in order to concisely refer to any parameter entering Eq.(1.2| 



3 



case, symmetry-violation is a relic from higher-dimensional operators at the non-symmetric vacuum, due 
to the truncation of the potential to dimension < 4 terms. To be definite, if at high energy, beyond 
a certain scale A, the symmetry holds, the potential V is then well approximated by its classical form 
(the symmetry-violating effects being negligible) and the Xf' at the scale A may be chosen as boundary 



conditions for the general parameters of Eq.(1.2), 



a; 



cl. 



A(A); 



, qm 







:i.3) 



such that V = V(A^^'(A)). At scales /i ^ A, however, symmetry-violating effects are no longer negligible 
so that non-trivial values of A^™ are generated by the renormalization group equations. 



We shall now enumerate possible symmetries one can impose to the potential of Eq. (1.2): 



Discrete Z^-symmetries: they are characterized by the transformations <I> i— )• e~*^*<^, where 
<I> = S,Hu,Hd and Qs,Hu,Hd the charges under the discrete symmetry group. They allow for 
significant selectivity among the complex terms of the general potential, while avoiding the prob- 
lem of an axion (unless the potential they induce is also accidentally {7(l)-invariant). Spontaneous 
breakdown of these symmetries (through Higgs v.e.v.'s) however generically leads to cosmological 
difficulties, in the form of a domain- wall problem pS) , which should then be addressed separately. 

1 . The complex doublet-terms are controlled by Qhu + Qh^ ■ Qhu + Qh^ = [n] causes no con- 
straint; for even n, Qhu + Qh^ = f [^] allows only for A5; other choices forbid all the corre- 
sponding terms. 

2. Complex mixing-terms are governed by both Qhu ~^QHd Qs- Only in the case {Qhu ~^QHd = 
0[n],Qs = 0[n]} are they all allowed by the Z„-symmetry. Otherwise, the relative choice 
of Qs and Qh,, + Qh^ constrains them, with the specific values Qs = ±{Qhu + Qtld) M) 
2Qs = ±{Qhu + Qh^) M and up to the exclusion of all these terms. 

3. The complex singlet-terms are governed by Qs, ranging from conservation of all {Qs = [n]) 
to exclusion of all, with the special cases 2Qs = 0[n], 3Qs = [n] and 4:Qs = [n]. 

A typical example for such a discrete symmetry and deserving particular attention is that of the Z3- 
symmetry with charges Qs,Hu,Hd — 1- this corresponds to the case of an underlying NMSSM.Invariance 
under $ 1— t- e 3 $ reduces the potential to the form: 



+ ml\S\^ + K^\S\^ + 



S'' + h.c. 



+ X'PlSflHul^ + XUS\^\Hd\^ + [AudSHu ■ Hd + X^'S*^H^ ■ + h.c] 
The tree-level conditions resulting from the NMSSM read: 



Ai = 

Ap 



9^ + 9" 



Xk 



As 



9 - 9 
4 

= kAi^ 



12 



; A4 — 
) A^i^d 



2 



A^A 



A2 
,2 



Ap , 



:i.4) 



:i.5) 



Our notations for the SUSY parameters follow those of 18 , except for the electroweak gauge 



couplings which we denote as g' and g for, respectively, the hypercharge U{1)y and the SU{2)l 
symmetry. 

Continuous global symmetries: here we mean essentially global phase transformations $ 1— )• 6**^*"$, 
that is C/(1)-Peccei-Quinn (P.Q.) symmetries 39 . Such symmetries are spontaneously broken by 



the v.e.v.'s of the Higgs fields so that they produce massless axions. They are also chiral in nature, 
so that anomalies will be generated at the quantum level (unless the field-content is enlarged so as 
to cancel them). Such symmetries are thus likely to stand only as approximate limiting cases. 
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1- {Qhu + QHd — 0' Qs = 0} is automatically satisfied: this is the hypercharge. 

2- {Qhu + Qh^ — ^jQs 7^ 0} preserves the doublet potential while constraining drastically the 
singlet couplings: 



V 



s-s 

PQ 



2 - (m?2^„ • Hd + h.c.) + ^\Hd\^ + ^\Hu\^ + Xs\Hu\^\Hd\^ 



A5 



Hu ■ Hctf' + (Ael-f^Mp + \y\Hii\'^)Hu ■ + h.c. 



+ + K^\S\^ + \l\S\''\Hu\^ + \i\S\''\Hd\^ + [Xm\S\^Hu ■ + h.c) 



:i.6) 



3- {Qhu + QHa 7^ 0) Qs = 0} constrains severely the doublet sector, as well as the mixing terms, 
while leaving the pure-singlet potential untouched: 



S-D 
PQ 



+ m||5|2 + K2|5|^ + 



+ + ^^-S' + ^>5|5|2 + + ^S^\S\' + h.c. 



2 3 3 ' ' 4 

+ A|,|5p|^„|2 + A^|5|2|i?rf|2 + [a^52|^„|2 + A^52|i?rf|2 + h.c 



:i.7) 



4- {Qhu + Q-ffd / 0,(^5 = -{Qhu + Q/^d)} is the 'usual' Peccei-Quinn symmetry (e.g. |40]) and, 
without loss of generality, one may choose [Qh^ = 1 = Qh^iQs = ~2). It induces a potential 
of the same form as that of the Za-symmetry (Eq. (1.4)), with the further requirement that As 



and Xp vanish. 



+ m||S|2 + k2|5|4 + A^|5|2|i/,|2 + xUS\''\Hd\^ + [AudSHu ■ + h.c] 



VpQ = m^J/7„|2 + m^J/7rf|2 + ^|Frf|^ + ^|i/j4 + A3|//„P|i/dP + A4|F„-Fd|2 



5- {Qhu + Q-ffd 7^ 0, = Qhu + Qh^j} is equivalent to the preceding case with the replacement 
S = S*. 

6. {Qhu + Q-ffd 7^ 0, Qs = \{Qhu + Q-ffd)} is a variant, concerning the singlet-doublet mixing- 
sector. This is again a subcase of the Za-potential (Eq, 1.4), with vanishing As and A^d'- in a 
coarse understanding of the term, this may be considered as the 'R-symmetric' potential. 



Ai, 



A, 



PQ, = mijjHuV + mij^\HdV + ^\Hdr + ^\Hur + X^\H^\'\HdV + XA\H^-HdY 

+ m||5|2 + K^\S\^ + A^|5p|i7„|2 + A^|5p|/7rf|2 + [a|/5*2f„ • + h.c] (1.9) 

Note that if one is interested in a SUSY-inspired model, this PQ'-symmetry would a priori forbid 
the XSHu ■ Hd term, resulting in vanishing tree-level conditions for most of the parameters of 
Eq.(1.9): it is therefore best understood as a R-symmetry at the SUSY level. 

7- {Qhu + Qh^ 7^ 0,(5s = —\{Qhu + Qh^)} is equivalent to the preceding choice, with the 
replacement S ^ S = S* . 

8. {Qhu + / 0, Qs / ± {O, 5, 1} {Qhu + Qho)] forbids all the complex terms, hence leading 
to another, more-constrained subcase of the Za-potential: 



V|q^ = m]jjHu\'' + mlj\Hd\'' + ^\Hd\^ + ^\HX + HHu?\Hd? + Xi\Hu-Hd\'' 

+ ml\S\^ + K^\S\' + Xl\S\''\Hu? + Xi\S\^\Hd\^ (1.10) 
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In the following, we shall focus only on Vpq and Vpg,, which both can be viewed as subcases of 

C/(l)'-gauge symmetries: they can be regarded as the gauged-version of the P.Q. symmetries, with 
the important consequence that the P.Q.-axion is now unphysical. They emerge naturally from 
C/(1)'-SUSY models, containing SM-singlets charged under the additional C/(l)'-gauge symmetry 
and breaking it spontaneously while acquiring v.e.v.'s. The simplest version of such models, with 
only one singlet, is called UMSSM |22| and leads back to the Zs-invariant Higgs potential, but with 
vanishing As and A|f, i.e. V^mssm ~ ^PQ- Eq. |l.8|. The further tree-level conditions are 



shifted from Eq. (1.5) according to (with Qs,Hu,Ha the Higgs charges under the C/(l)'-symmetry and 



gz', the coupling constant): 

Note that the SM-fermion sector is also charged under the U (l)'-gauge group, so as to ensure 
invariance of the usual Yukawa terms. To avoid a chiral anomaly of the 11(1)' symmetry, an exotic 
fermion sector will also be necessary. 

One may also write tree-level conditions of a different form, not protected by any symmetry: this is 



for instance the case in the nMSSM, where a or a Zi^ symmetry |2l| is imposed at the level of the 
superpotential, so as to forbid all renormalizable pure singlet-terms, then broken explicitly by gravity 
effects, in order to arrange for an effective tadpole term (so as to break the resulting P.Q. symmetry), 
broken also explicitly by the soft-terms. The tree-level potential then differs from the Z3 case ( |1.4[ ) by 
the requirements: 

X^^ = K = As = ; At, / (1.12) 

We hence define: 

= mjjjH^ + mj^jHdl^ - {ml^Hu ■ Ha + h.c.) 

+ + [XtS + h.c] + A^|5|2|F„|2 + Af>|5p|Fd|2 + [AudSH^ ■ Ha + h.c] (1.13) 

While the absence of a residual symmetry at low-energy is a deliberate feature of the nMSSM (in order 
to circumvent both axion and domain- wall problems), the resulting lack of protection of the tree-level 
couplings at low-energy will lead to sizeable consequences for the parameter reconstruction at the loop- 
level, as we will see later. 

1.3 Mass matrices 

Spontaneous symmetry breaking is achieved when the scalar fields develop a v.e.v., 

{Hu) = f ° ) , {Hd)=( ), {S) = s (1.14) 



Imposing the minimization conditions associated with the most general potential in Eq. (1.2), one may 
trade the parameters rrij^^, n^'jj^, rrig for the v.e.v.'s Vu, Vd, s. Introducing the usual definitions v = 



6 



' v'^ + '^J — 174 GeV, tan /3 = Vu/vd, we can write these relations as 



ml 



And + Aud + {\p + \p + Xm)s 



2Ads + {\% + 2~4,)s 



XlV Cp - (A3 + A4 + A5) V + [\&V Sp - 7^12) tp + 3w S2/3A 



Aud + Aud + {\¥ + + \m)s 



St: 



2Aus + (A^ + 2A^)s 



A2'V^s| - (A3 + A4 + A5) t>^c| + (Atv^c^ - mfs) ^/3^ + 3?;^S2^A6 



(1.15) 



And + Aud + 2(A:^^ + + AAf )s 

(A,, + 2A^s)4 + (Arf, + 2Af,s)ci 



As + As + {2k^ + 4 + '^l)* 



2s 

\m22 \d22 '^T 2 
XpV Sp — XpV Cp /i^ 



The quadratic terms in ^ provide us with the charged Higgs mass matrix: 



M 



[Aud + Aud + (Ap + \p + \m)s)s 



;(A4 + Xb)s2P - - XjCp V - 



(1.16) 



Its diagonalization expectedly delivers (massless) charged Goldstone bosons = cos/3 — sin/3 ffj^ 
and the physical charged Higgs = cos/3 -fT^ + sinfiH^ , with mass: 



S2/3 



{Aud + Aud + {X¥ + + Xm)s)s 



12 



We turn to the CP-odd squared mass matrix, written in the basis (aJJ, a„, a^): 



(A4 + A5)s2/3 - XqsI - A7c| ) - 
„o ^0 



(1.17) 



MI22 



Aud + Aud + (Af + A|f + AAf)s j s + (Ags^ + Ayc^ " ^552/?) 
Aud + + (Ap + \p + AAf)s) s + {\sl + A74 - A5S2/3) 
A«d + + 4(A|/ + A^^)s 



V — m 



12 



2 2 



l'^S2/3 



2s 



A, 



3^5 + ^ + (44 + kI)s 



(1.18) 



V 



--24-- 



M 



P13 



Aus + 4A^sj + (^Ad, + 4A?,sj c; 
A«d + ^d + {X^P + ~Xp + Aa/ )sj s + (A6s| + A74 - A5S2/3) 
Aud — A„(i — 2(Ap — Ap )s 



2 2 



P23 



A 



ud 



Aud - 2(A*^ - Af )s 



The neutral Goldstone boson = cos/3aJ] — sin/3a^ can be isolated through the rotation with angle /5 
and we are left with the 2x2 matrix A^p, in the basis (a^,, a^), with a^, = cos/3a^ + sin/3a^ 

2 



P' 11 



S2/3 



A^.d + A„rf + (A^^ + Af + AAf)s ) s - (A5S2/3 - Ae^ - Ayc^) 



^2\„,2 



V — m 



12 



P' 22 



Aud + A„rf + 4(Af + A|f )s 



f^g2/3 

2s 



3As + ^ + (44 + K^)s 



(1.19) 



P' 12 



+ 4A^sj sj + l^Ads + 4Af,sj 4 
A-ud — Amc( — 2{Xp — X^{)s 



V 



--24-- 



^We use the shorthand notations C/s = cos/3, S/3 = sin/3, S2,9 = sin2/J, = tan/3 etc. 
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A^p/ is diagonalized in the subblock of the physical states {a^,a^) by the orthogonal matrix P' , to give 
the two physical CP-odd squared mass m^oj ^'^oi such that 



^1 '^2 
,2 2 



diag{mlo,mlo) = P'Ml,P' 
Finally, the CP-even squared mass matrix, in the basis {h^, h^, h^), reads: 



(1.20) 



-Mill 

-^522 

-Ml 33 
-M|i2 

^Il3 
^|23 



Aud + A^d + {\p + \p + \m)s ] s + (Xesj - SXjcj) 



:m 



^2\„,2 



V — m 



12 



Aud + Aud + (A|f + A;^ + Xm)s ) s + {Xjcj - SXqs 



[M 



,2\„,2 



V — m 



12) 



tn + 2Xiv'^cl 



tl^ + 2X2v'^sl 



Aud + A. 



ud 



2s 



+ 



As + As + 2(2^2 + 4 + 4)s s - {Auss} + Ads4 



Xt 



And + Aud + {X^P + ~Xp + Xm)s s + [(A3 + A4 + A5)s2/3 - 3(A64 + hcfj] + m?2 



^nd + Ad + 2(Af + Af + AAf)s ?;s;3 + 2 A^, + (A^ + 2A^)s 



Ad + Ad + 2(A^^ + Af + Am)s vcfi^2 + (A^ + 2A^)s 



a.2i) 



which is diagonalized by a 3 x 3 orthogonal matrix S, resulting in three CP-even squared masses m} 



m^O) "i^O) such that 



diag(m^o, m^o, m^o) = SM\S ^ 



(1.22) 



We are thus finally left with seven physical Higgs particles, once the three Goldstone bosons G", G^, 
giving mass to the and Z'^ bosons, have been discarded. In the particular case of the C/(l)'-gauge 
symmetry, however, the P.Q.-axion (associated to the vanishing eigenvalue of A^p/) is also unphysical 
(giving mass to the Z'-boson, gauge-field of the J7(l)' symmetry |22|), so that we are left with only one 
CP-odd physical mass. 



2 Reconstruction of the effective parameters 
2.1 Masses and mixing angles as physical input 

From an experimental point of view, the 'Aj' parameters are not directly accessible: they will enter as 
combinations within the expressions for the Higgs masses and self-couplings. The latter can hopefully 
be accessed through the experimental measurement of physical quantities. 'Inverting the system', we 
can therefore trade some Aj parameters for such physical input. In the simplest case, one would directly 
use the Higgs masses and their mixing angles, assuming these can be measured (e.g. from fermion/gauge 
couplings), as the new, physical input. For the 2-doublet-|-l-singlet system, such quantities provide us 
with 12 conditions (input measurements) on the Aj's: the masses of the 2 CP-odd bosons, 3 CP-even 
and 1 (complex) charged Higgs; the mixing angles from the CP-even (3), CP-odd (1) and the Goldstone 

(1: /3) sectors; finally, the electroweak v.e.v. v = \Jv'^ + v'^ (from M\y for example). Should those 
twelve relations prove insufficient to determine all the Aj's (as is obviously the case for the most general 
potential), one would have to resort to Higgs self-couplings (or input from another sector) in order to 
fully determine the parameters. Accessing such self-couplings would require that double or triple Higgs 
production are kinematically open. This task would most comprehensibly done at future linear colliders. 
In the meanwhile, the measurements of masses and mixing angles still allow for a partial inversion. 
We will assume in the following that all the Higgs-masses have been measured. Note that this hypothesis 
is somewhat optimistic since singlet-like fields do not couple directly to SM-fermions and gauge-bosons, 
hence are essentially elusive: only when there is substantial mixing with the doublet states can we expect 
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to access them without having to rely on multi-Higgs couphngs. As for the mixing angles, assuming all 
the Higgs states have been observed in SM decay-channels, one can derive them from the couplings to 
fermions (note that leptonic decay channels are likely to give cleaner information) and gauge bosons. For 
a type II model, we have (taken from [18]): 



V2 
Y 



<^%t% = -^^cpPl, (2.1) 



V2 



Y 

H+hL^k = Ytcp 
H-tLb'k = -Ybsp 

VSp VCfS vcp 

and (we mention here only the 1-Higgs to 2-gauge couplings; note that, albeit more difficult to measure, 
2-Higgs to 1-gauge as well as quartic couplings shall play a very important role for testing the model): 



where. 



V2 



i2 



n2 



hW+W- = gf,u^v{cfsSa + sisS,2) (2.3) 

Combining Higgs couplings to the vector bosons with those to up/down fermions, one can access e.g. 
Sii/Si2- Moreover, one may be tempted to use Higgs decays into two photons to extract information 
about the mixing angles: even admitting that such processes are dominated by quark loops, the corre- 
sponding relation of branching ratios to mixing angles is already non-trivial and would require an involved 
extraction procedure for exploitation. 

Unitarity relations could also prove useful. For example, a 'missing' matrix element Sij could be recon- 
structed from 

3 3 
""^SikSjk = 6ij = ''^SkiSkj i,j = 1,2,3 (2.4) 

fc=0 k=0 

A possible (naive) strategy to reconstruct the mixing angles would be the following: having measured 
the charged Higgs decay into third generation quarks, one could then deduce t^, since the ratio mtfi/v is 
fixed by SM measurements. Then the (doublet) elements Sji, Sj2, P'n could be obtained unambigously 
from the decays of neutral higgs states into fermions and gauge-bosons. The unitarity relations would 
finally provide the magnitude of the and P/2 (singlet) elements. 

Note finally that, while a full experimental determination of the Higgs mass matrices may seem over- 
optimistic in the short run, there exists a practical case where we have access to such data: it is that of 
the output of spectrum generators (e.g. the publicly available NMSSMTools, |31j). We will resort to that 
practical application in the last part of the present paper. 
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2.2 Partial reconstruction in the general case 



Considering the general potential of Eq. (1.2) and discarding any assumption as to an underlying model, 
a complete reconstruction of the 29 parameters (28 of which being relevant) cannot succeed with only the 
twelve mass/mixing conditions, hence calls for the measurement of Higgs self-couplings. Yet, information 
from Eqs.( [rT7ll.20|1.22D can already be implemented in a partial reconstruction: 



([{A 



ud + Xp s)s 



m 



12J 



S20 



-3^5. + {A^d + 4A^s) ^ - ^ = ml,P:^ + Xl 
^[iA^d + X^s)s->fv^S2(s 
[{A^d + X^h)s ■ 
[{Aud + Xj^s)s ■ 



m 



12 



m 



12 



m 



\tp + 2X,v\l = ml,Sl-Xl 
\t-^ +2X2V^ si = mloSj. 



Ass + 4k2s2 + Aud^,S2p 



A| 



\2 



-{Aud + X^Js)s + (As + X4)v^S2/3 + = mloSiiSi2 - Xf 



-{A^d + 2Xfs)vsp + 2X%svcp = mf^oSnSis 



Xl3 



-{Aud + 2Af s)wc^ + IX'psvsf^ = mloSi2Si3 - Xf 



(2.5) 



1 \ 1,2,3 \12 \ \ 1,2,3 ,12,13,23 • i 

where Ap , Ap , A± , A^. , A^. are given by 



AF 

XI- 



\23 



{And 



2 ( + (A^j + Aa/ )s ) s - (A5S2/3 - Aes 
p 

M 



2XIjs)v 



And + 4A 



2s 



+ 



t + (44 + 4)^ 



s + 



A7C2)1 



{Aus + ^Xls)sl + (^d, + 4A^s)c2 



+ 24 



2 

S2/3 



(i„d + (A|f + Aa/)s)s - IA5S2/3 



Ae^ 



A74 



+ {Xf + Xm)s s + (A6s| - 3A7c2 )w2 



+ (A^^ + Aa/)s J s + (A74 - 3A64)^^^ 

Aud 2s^ ^•S'S + 2(k|. + /i|.)s^ — ^A^sS^ + AdsCjj 

-{Aud + (X^p + Xm)s)s + A5S2/3 - 3(A64 + A74 

vsp + 2(2Af,s + Ads)vcp 



Aud + 2(A|f + Xm)s 
Aud + 2(A|f + Aa./)s 



vcp + 2(2A^s + Aus)vsii 



(2.6) 

Our (arbitrary) choice in ordering the parameters within Eqs.( 2.5|2.6 ) was guided by the terms that are 
relevant at leading order in the n/NMSSM and the UMSSM potentials: beyond Tn?u^ 
are common to the three models, those are given by 



m 



Hu' '"-5' 



?r, which 



NMSSM 
nMSSM 
UMSSM 



2 ^"''^ At, a 



Al-4 

Ai_4, m^2-> X 

u,d 



d ' 
Al-4, Ap , K^, Aud 



ud, As 
ud 



These parameters were collected on the left-hand side of Eq.(2.5), while the remaining ones enter the 
right-hand side through Eq. (2.6). 



Note that the relations of Eq.(2.5) hold at any order (since Eq.(1.2) is the most general renormalizable 
potential satisfying the gauge-symmetry). A practical use of Eq. (2.5) would lie in a model-independent 
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analysis of a 2-doublet+lsinglet potential (in order to discriminate among models, constrain them through 
precision tests). Then the twelve mass conditions can be used to simplify twelve (arbitrarily chosen) 
parameters, hence leaving the remaining couplings as the relevant degrees of freedom intervening in / to 
be determined from the Higgs self-couplings. Not much predictivity should be expected, however, in this 
general case. 



2.3 Reconstruction at the classical level in the constrained models 



We focus here on the specific cases inspired by the SUSY models: V-f, Vpq, Vpq, and Vf^. Note that 
such potentials are considered at the classical order: quantum effects and explicit /spontaneous breaking 
of the symmetries in principle destabilize those potentials to generate the most general one. At this 



leading order, however, the Eqs.(2.6) vanish, leaving Eqs.(2.5) in a very simple form. Note additionally 



the further requirements for each potential: 





: m?2 = 


= Xt 


= 




: As = 


K^ = 


- \P — 







■■ As = 




= m\^ 


= At = 


•^PQ' 


■■ As = 


A-ud 


= "li2 


= At = 



We end up with eleven classical parameters and eleven condition^for both the potentials Vpq and Vpqi. 
In these cases, all the parameters in the Higgs potential can thus be reconstructed (at leading order) 
from Higgs masses and mixings: this procedure is explicitly carried out in appendix [A] Eqs.( A.5|A.6 ). 
In the case of Vf^^ , the thirteen classical parameters cannot be fully determined from the twelve conditions. 
The remaining degree of freedom is conveniently chosen as the singlet v.e.v. s: the reconstruction is 
also given in appendix |Aj Eqs.( A.l|A.2 ). Several tracks can be followed in order to determine this 
remaining degree of freedom. The first one, sticking to the Higgs potential, would consist in relying on 
trilinear couplings, such as h^H^ H~ or h^a^a^, where the neutral Higgs fields would be largely singlet in 
nature: kinematical limits and the elusive nature of singlets would tend to disfavor this strategy. Another 
possibility would be to input information from some other sector (if any): measurement of the higgsino 
masses in the NMSSM could provide the missing information. Finally, a more predictive option would be 
to enforce some additional requirement, such as relations among the tree- level couplings: the tree-level 



relations of the NMSSM, 



A" 



1 or 



{AM)2-(a+6) 



= 1 (where a, 5 are real numbers), for instance, or a 



aX^+bXj, A 



J, may be used as guidelines. 



measure of the P.Q. symmetry breaking, such as 

Finally for V^, we have twelve parameters and twelve conditions. Yet a full inversion is not possible 
either, because CP-even and CP-odd singlet masses are explicitly degenerate in this potential, leaving a 
bound system. The remaining degree of freedom is again chosen as the singlet v.e.v. s in appendix lAl 



Eq.(A.7), but could be replaced by e.g. At, as a measure of the violation of Z3, for example. 
So far, we have considered only the Higgs potentials separately. Moving explicitly to the underlying 
SUSY models, however, the Aj's are further constrained by the tree-level relations resulting from their 
supersymmetric origins: we count 7 parameters in the nMSSM Higgs sector (At, ?7ii2, 



m. 



A, A),), 7 in the NMSSM as well (ml 



A, Ax, K, Af,) and 6 in the UMSSM {mjj^, 



m 



m 



5' 



A, Ax, gz'', note that we regard the Higgs charges under U{1)' as fixed). Those parameters are then 



over-constrained by Eq. (2.5) and one should thus consider the remaining conditions as a measurement 
of the deviation from tree-level conditions due to higher orders (we remind here that the tree-level 
relations induced by the model of origin among the parameters of the potential are likely to be spoilt 
by quantum corrections). Depending on the information at our disposal in the remaining spectrum (e.g. 
SUSY masses), such conditions may be used to estimate the missing parameters (e.g. sfermion masses or 
trilinear soft couplings) or regarded as precision tests of the model. Note that if the SUSY spectrum is 



'The explicit presence of a P.Q.-axion, identified as a?, leads to one trivial condition in the CP-odd sector. 
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sufficiently documented as well, this measurement of the Higgs parameters at leading order, would allow 
for a (perturbative) computation of all the Aj's within the specific models at higher orders. 



2.4 Reconstruction at the loop level: NMSSM vs. nMSSM 

Now we want to apply this formalism to higher order effects. The purpose is simple: it has been shown 
that, in the MSSM, the bulk of the corrections in Higgs-to-Higgs couplings could be absorbed in writing 
such couplings in terms of the corrected masses (see for example [4l||42| and the third reference in |43j ) ; 
could a similar recipe apply to the 2-doublet+lsinglet setup? A first strategy is the one presented at the 
end of the previous subsection: in a definite model, the Higgs spectrum may allow for a determination of 
the Higgs parameters at leading order; then, provided sufficient information from the other sectors stand 
at our disposal, reconstructing all the Aj's at higher order is simply a matter of perturbative calculations. 
Yet, this approach relies on a heavy machinery and on input which is external to the Higgs sector. We 
would like to consider cases where input from the Higgs sector only (or almost only) would already 
improve on the simple tree-level expression for the Higgs self-couplings. 

In principle, whatever the potential looked like at the classical level, quantum corrections will generate 



contributions to all the parameters in the general potential - Eq.(1.2) - (unless a symmetry protects cer- 



tain parameters, but we have seen that such symmetries are spontaneously broken by the Higgs v.e.v.'s 



anyway). Therefore, while the partial-inversion of the general case (Eqs.( 23|2.6 )) is still possible, little 
predictivity or practical use is to be expected from such relations, because the number of undetermined 
parameters is high. To extract meaningful information, beyond the leading order, from the Higgs spec- 
trum, one would need the corrected potential to retain a sufficiently simple form beyond the classical 
order. 

To be more specific, we consider a tree-level potential of the form {H representing any of the Higgs fields, 
/x^, a bilinear. A, a trilinear, and Aj, a quartic coupling): 

Vtree = fJ^'^H^+AH^ + XH^ (2.7) 

We now include the radiative corrections, which shift the potential as: 

Veff = ifi^ + 6fi^)H^ + {A + 6A)H^ + (A + 5\)H^ + S^^H"^ + 5AH^ + 5\H^ (2.8) 

where Sji'^, 5A and 5\ represent corrections to parameters existing at tree-level, while 5A and 5\ 
denote new couplings which were forbidden by symmetries at tree-level and emerge only at the radiative 
level. Neglecting numerical coefficients, the corrected Higgs mass and the trilinear self-coupling g will 
read (schematically): 

fm^ ~ + + 5 + {A + 5A + 6A){H) + (A + 5A + 5~X){H)^ = ml^^ + 0(6, S) 
[gc^A + 6A + 6A + {\ + 6\ + 6\){H) = gtree + 0{5, 6) 

(with the short-hand notation 6/6 for loop induced corrections to parameters present /absent at tree level.) 
We now assume that we have access to the mass m^, either from experimental data or from a spectrum 
generator. Using gtree in the computation of physical quantities (branching ratios, cross-sections) will 
result in an error of order O (^^). If we use the expression for the corrected mass to inverse (partially) 

the relation between mass and tree level parameters, we obtain: 6 = 5^2 -|- 0{6), where (5„2 symbolises 
the result of the inversion procedure. The trilinear couplings then provide: g^2 = g + C(<5)) resulting in 
an error of O at the level of cross sections/branching ratios. Claiming that the inversion procedure 
carries any improvement with respect to a simple tree-level evaluation holds at the sole condition that 
radiative corrections 6 to tree-level parameters are more important, in magnitude, than the contributions 
6 to other operators. Otherwise, even if we identify the parameters subject to large contributions, it is 
unlikely that the mass-matrices would suffice in determining both these parameters and those intervening 
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at tree-level, unless we input some additional tree-level relations, as in the case of the matching conditions 
in Eq.([L5]), dTTTj ) or ^J^. 

This discussion shows that, to extract some benefits - beyond the leading order - from the conditions 
relating masses to effective parameters, we need to identify which terms are potentially subject to large 
radiative corrections. A simple criterion can be invoked at the one-loop level: it is that of the leading 
logarithms. To identify those, we simply resort to the Coleman- Weinberg [28] one-loop effective potential 
and analyse the outcome for the special case of the SUSY-inspired models under scrutiny. This method 
has long been employed for the computation of corrections to the Higgs masses, both in the MSSM [43] 
and in the NMSSM T8|31[ (and references therein). In this approach, the effective corrections to the scalar 
potential at a scale A are determined by the field-dependent tree-level mass matrices M|(5, H^, Hu, ■ ■ ■) 
of the various fields $ entering the spectrum, according to (in the Z)i?-scheme, but note that we shall be 
interested in the logarithms only): 



647r2 



(2.10) 



Here C$, which counts the degrees of freedom, takes the values 1 for real scalar fields, 2 for complex ones, 
—2 for Majorana fermions, —4 for Dirac fermions and 3 for real gauge-fields. Note that we are interested 
in the Higgs potential solely, so that we will retain dependence on S, Hd, Hu only, within M^. Moreover, 
we consider no EW- violating effects so that we will not expand the doublet fields H^, Hu around their 
v.e.v.'s (except within logarithms). Additionally, the S'C/(2)£,-symmetry can then be invoked to retain 
only the neutral Higgs fields S, H^, Hu (the dependence on the charged Higgs fields can then be restored 
afterwards in virtue of SU{2) l: only the A3 and A4 parameters cannot be disentangled in this fashion, but 
both parameters being present at tree-level in the models we consider, this will be of little consequence 



for our analysis). We then determine the contributions to the parameters of Eq.( 1.2 ) by letting the singlet 
take its v.e.v., S = s + S, then truncating Eq.(2.10) to renormalizable terms, finally projecting on the 



couplings of Eq.(1.2). 



The results of our analysis of the large logarithms, in the cases of the NMSSM and nMSSM, are provided 
in appendix |Cj The situation of the NMSSM is quite simple: leading logarithms favor Zs-conserving 
terms. We can thus claim, for this model, that the inversion procedure for the Zs-conserving potential, 
presented in the previous subsection and Eqs.( A.l|A.2 ), improves on the tree-level implementation of the 
couplings and actually includes leading-logarithms automatically. Note that, as defined in Eqs.( A.l|A.2 ), 
the effective Zs-conserving parameters are directly determined in terms of physical quantities, meaning 
that they do not depend on the renormalization scale A : they are simply the parameters of the effective 
Zs-conserving potential associated with the physical Higgs spectrum. What we checked explicitly in the 
Coleman- Weinberg approach (which depends on the renormalization scale A) is that this constrained 
form of an effective potential was legitimate at least up to leading logarithms. Beyond, the effect of the 
Zs-violating terms (due to the truncation of the potential to operators of mass- dimension < 4) cannot 
be neglected. In the case of the nMSSM, however, potentially large logarithms affect non-classical terms. 
In fact, all the sectors contribute to the Za-conserving parameters (including those vanishing at tree 
level in this model). Additionally, logarithms originating from the nMSSM Higgs sector (the only sector 
which is sensitive to the breakdown of Z3 at tree-level) also affect Zs-violating terms. Inclusion of the 



leading higher-order effects from the inversion procedure of subsection 2.3 Eq.(A.7), thus seems dubious 
in this case. It seems natural to ascribe this difference of behavior, between nMSSM and NMSSM, to 
the protection of the parameters by the Za-symmetry, which albeit spontaneously broken by the singlet 
v.e.v., continues to favor Zs-conserving terms within the NMSSM. We should thus expect a similar 
property, whatever the Zs-symmetric model is (SUSY or not), and, beyond the Za-symmetry, in any 
model retaining a symmetry (or approximate symmetry) at low-energy, e.g. PQ or PQ' . 
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3 Phenomenological consequences for the NMSSM 



We now explain how, with the formalism derived above, we can improve the computation of some ob- 
servables in the NMSSM. As we have already highlighted, there is one practical case where the Higgs 
spectrum is fully available: it is that of the spectrum generators. The Higgs masses are often, in such 
a case, corrected, while couplings are typically taken at tree-level. In the case of the NMSSM, we have 
shown that leading quantum corrections could be absorbed within the tree-level parameters of the Higgs 
potential. This allows us, at a very cheap cost, to improve the accuracy of the Higgs self-couplings by 
reexpressing them in terms of the Higgs masses and mixing angles provided by the spectrum generator. 
We refer to appendix [B] for the explicit expressions of the couplings in terms of the effective parame- 



ters. We implement this recipe both within NMSSMTools_3 . 2 . 31 32 and within SloopS 35 36 , and 
investigate phenomenological consequences. 

3.1 Impact on Higgs-constraints within NMSSMTools_3 . 2 . 

NMSSMTools_3 . 2 . includes several phenomenological constraints on the NMSSM parameter space, orig- 
inating e.g. from LEP [44j , TeVatron |45j , B- and T-physics |46] as well as early (now outdated but in 



the process of getting updated) LHC-data 47 . The Higgs-sector evidently plays a central part in the 
interplay of these experimental limits and we would like to investigate whether our analysis could have 
meaningful consequences at this level. 

The basic routine mhiggs . f of the NMSSMTools Package computes the corrections to the Higgs mass 
matrices, incorporating typically leading-logarithmic effects (although leading two-loop contributions 
from the fermion sector are also implemented) at a scale determined by the stops and sbottoms, then 
rescaling the fields at the EW scale, finally adding pole corrections (this whole procedure is more precisely 
described in the appendix C.3 of |18j). In this subsection, we shall be relying on this routine for the 
calculation of the Higgs masses and mixing angles from the NMSSM parameter input. Note however that 
NMSSMTools offers a second possibility, which consists in the evaluation of the Higgs masses according 
to j48j , including the full one-loop corrections as well as the two loop 0{atas + 0ib0.s) (with at^b = ^tb/^'^) 
effects in the effective potential approach. This option will be used in the next subsections. Whatever 
the source of the masses and mixing matrices however, we will treat the latter as input for the physical 
Higgs matrices, allowing us to compute the A^'s. 

The Higgs couplings implemented within NMSSMTools_3 . 2 . are actually not purely tree-level couplings: 
possibly large radiative corrections from the quarks of the third generation are included as well, as 



explained in the last paragraph of the appendix A. 2 of 18 . One can check that these corrections arise 
from (s)fermion contributions to Ai^2 (as one can recover considering our results for the Coleman- Weinberg 
analysis in appendix [C|): such effects are thus in principle automatically incorporated within our procedure 
(given that the corresponding contributions to the Higgs mass are also implemented within NMSSMTools). 
The choice of the singlet v.e.v. s deserves an additional comment, since it cannot be extracted from 
the masses. After comparing with a few variants leading to minor deviations (a few percent) at the 
numerical level, we settled for the simple definition s = ^es/^i with //gfr and A the parameters inputed 
in NMSSMTools. Note that this choice is coherent with the recurrent use of /ieff/A as the singlet v.e.v. 
within the routines of NMSSMTools. 



To perform the comparison, we simply implement our corrected Aj's (see Eqs.( |A.l|A.2| )) and the ensuing 



triple Higgs couplings (Eq. ( B.7lB.8|B.9" )) within the routine decay . f of NMSSMTools, computing the Higgs 



decays. A flag enables us to choose between the original setup of NMSSMTools and our modified version, 
which we denote in the following as NMSSMTools*. 

Admittedly, the modification is essentially a fringe effect and one needs to go to a region of the parameter 
space where the Higgs self-couplings intervene very finely to discover substantial deviation between the 
two approaches. We thus consider a specific region in the NMSSM parameter space, characterized by a 
light CP-even Higgs with mass typically under 100 GeV, sizeable singlet-doublet mixing Sf^ ~ — 100% 
and a light CP-odd Higgs with mass < 10.5 GeV, allowing for — t- 2a5 decays. Such a scenario is 
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Figure 1: Constraints in the planes defined by (m^o,Sf^) (left-hand side) and {mf^o,Sf^) (right-hand 
side) for /ies = 300 GeV. Black dots correspond to the points on which we perform the scan (without 
collider constraints); yellow dots are allowed by NMSSMTools_3.2.0 while red dots signal points allowed 
with self-Higgs couplings defined as in our procedure NMSSMTools*. 



possible e.g. in an approximate Peccei-Quinn limit ^ ^ 1, with the parameters (NMSSMTools input: refer 
to jsT 32 ; we use the index 'sferm' to denote any of the sfermions): tan/3 = 5, A = 0.5, k = 0.05, 

Aferm., fJ'es £ [100; 900] GeV, \A^\ < 30 GeV and Ma S 

1 or 2, depending on the specific 



6Mi = 3M2 = Ms = 1.2 TeV = mrferm. = 
[0; 4] TeV. Incidentally for those points, the doublet-like state {i 



point) reaches a mass of ~ 125 GeV in the limit of singlet-doublet decoupling ~ 0. Note however that 
we did not specifically attempt to preserve this feature of a Higgs state at 125 GeV (so that for significant 
mixing, we may have typically m^o ~ 90 GeV while m^o ~ 150 GeV): we simply mean to show that 
our procedure is liable to affect the output of NMSSMTools. Note finally that, for simplicity, we discard 
constraints from {g — 2)^, which may be taken care of separately, by tuning the slepton sector. All other 



collider constraints 
are kept. 



from LEP, i3/T-physics, TeVatron or early LHC data - within NMSSMTools_3 . 2 . 




III hi ( GcV) 




1,1 (OcV) 



Figure 2: Constraints in the plane defined by (m;j0,5^3) for G [100; 900] GeV. On the left-hand side, 
the scan with the yellow dots uses tree-level Higgs couplings, while the corresponding one on the right-end 
side is obtained with the couplings implemented in NMSSMTools, adding fermion corrections. The color 
code is otherwise similar to that of the previous figure. 

We first specialize to the case /ies = 300 GeV: the lightest CP-even Higgs is then dominantly singlet 
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(5^3 ~ 70 — 100%). We display our results for this scenario in FigjT] The black dots represent points 
on which we scanned (with no collider constraints applied; note that their distribution is an artifact of 
the scan and should not be paid particular attention), the yellow dots, the output of NMSSMTools with 
the original Higgs couplings. Our results (NMSSMTools*), with corrected Higgs couplings, are depicted 
by the red dots. Since we scan over two variables, the output is two-dimensional in the (m^o, 5'^3)-plane. 
In the plane (771^0,5^3), the constraint m^o < 10.5 GeV <C mj^o reduces the apparent dimensionality to 
1. To investigate the whole {nif^o, 5'^3)-plane, one may additionally scan on the parameter fi^g, which we 
show on Figj2] There, however, the plot on the left-hand side corresponds to the case of purely tree-level 
couplings (for the yellow dots), which we obtained by removing the fermion-corrections from the original 
couplings implemented in NMSSMTools. The obvious conclusion is that, although partially compatible, 
the yellow and red dots do not exactly coincide, so that the details of the constraints are affected by our 
procedure. Note that both points admitted by NMSSMTools while excluded by NMSSMTools* and points 
admitted by NMSSMTools* while excluded by NMSSMTools are to be found. 

We insist however on the fact that such displacement effects in the acceptable points of the parameter 
space are noticeable only because we considered a region where phenomenological constraints on the 
Higgs spectrum and decays are particularly severe, due to the presence of very light Higgs-states which 
need to be sufficiently 'invisible' to escape experimental limits. A slight perturbation of the Aj's is then 
liable to result into insufficient — )• 2a5 ('invisible') branching ratio, excessive decays into e.g. SM- 
fermions and/or excessive a^-signals in S-physics: in such cases, an increased 'invisibility' of the light 
Higgs states, that is increased singlet components, is required. If, on the contrary, the perturbation of the 
Aj's stabilizes the invisibility of the light Higgs states, additional parameter space becomes available. Given 
the sensitivity of such regions to perturbations and the complex interplay of constraints at stake there, 
it is very difficult to predict to which extent the allowed parameter space would be shifted or not. In 
any case, a detailed new analysis on the NMSSM parameter space is beyond the scope of this work. 
Our procedure simply ensures the consistency of the calculation at leading-log order, the couplings being 
adequately related to the spectrum. 



3.2 Implementation in SloopS 



In SloopS, the complete spectrum and set of vertices are generated at tree-level from the NMSSM 
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There, g, g 



Vd are determined by 



SUSY and soft - parameters through the LanHEP package 
the physical input Mz, Mw, v and tan/3. Then the radiative part of the Higgs potential needs to be 
implemented: the tree-level Higgs parameters A^, given in Eq.(1.5), are thus shifted as A, = A^ -|- AA,, 



AAj defining the radiative corrections to the parameters Aj of Eq.(1.4). Yet, the corrected Higgs masses 



are not computed within SloopS, but imported from NMSSMTools through the SLHA interface. Applying 
the inversion procedure (Eqs.( [A!l|A.2[ )), we obtain the Za-invariant AAj's from the inputed masses, 
diagonalizing angles and tree level A^'s. From now on we will call this procedure the 'effective physical 
potential approach' and refer to it through the acronym 'PhA'. The complete set of Feynman rules is 



then derived automatically in the FormCalc 50 conventions, the latter performing the calculation of the 
decay width. 

A powerful feature of SloopS is the ability to check gauge invariance of results through a generalized 
non-linear gauge fixing, which was adapted to the NMSSM |37j. The gauge- fixing Lagrangian can be 
written in a general form 

Cgf = -^F+F~ -^{F^l^-^lF^'l^ (3.1) 
c,w '^iz 24a 
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where the non-hnear functions of the fields F are given by 



F+ = {d^-ieaA^-igcwpZ^W^' + 

+iiw^ {v + 5iHi + 52H2 + ^3^^3 + i{f^G^ + PiM + hM)^ 



(3.2) 



F^ = d^A^" 



The parameters a, /3, - ■ ■ , €3 are generahzed gauge-fixing parameters. We also set S,a,z,w = 1 to keep a 
simple form for the gauge boson propagators. 

The ghost Lagrangian Cch is established by requiring that the full effective Lagrangian is invariant under 
BRST transformations. This implies that the full quantum Lagrangian, with Cc the classical Lagrangian, 

= ■^C + ^GF + C-Gh (3.3) 

be such that (5brs>Cq = and hence 5brs^gf = —^^^is^Gh [51j . The BRST transformation for the gauge 
fields can be found for example in [Sl] . The NMSSM specific transformations for the scalar fields can be 
. For the decay — )• 77 not all the parameters are relevant: only a and 5i are. 



found in 
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3.3 The decay 77 

The diphoton decay is an interesting process to investigate due to the relevance of this channel in the 
recent discovery at LHC and because the gauge invariance is fully at play there. Indeed, in the SM, 
the W-boson loop, together with the top-quark one (the latter being of course gauge invariant, as the 
remaining contributions), dominate the decay. The calculation of the diphoton rate in the non-linear 



gauge was originally performed in 52 in order to simplify the calculation of the Higgs decay into two 
photons in the SM. In short details, with the specific choice a = —1, the W^G^^ coupling vanishes and 
the gauge-boson loop calculation is made easier. In our calculation we will refrain from adopting this 
choice in order to preserve the ability of checking the cancellation of the unphysical gauge-dependent 
part in the gauge loops. We will discuss this effect only at one-loop order, meaning that we consider 
only the LO decay width: this will be sufficient for our purposes since we do not aim at a more precise 
evaluation. Nevertheless it is worth reminding that the full two-loop EW+QCD corrections for the SM- 
like Higgs decay into 77 are known and under 2% [53] below the WW threshold. The full two-loop SUSY 
corrections are as yet unknown and anyway, our procedure would not be suited for such a calculation as 
the renormalization of the tree level Higgs potential would be mandatory whereas our renormalization is 
effective and explicitly breaks SUSY. We only aim at showing how one may consistently use the radiatively 
corrected Higgs masses for an improved LO calculation of this decay. 

Once we trade the "Aj" parameters for the masses, using Eqs.( A.1|A.2 ), we can re-express the Higgs 



self-couplings, obtained from the restricted Zs-invariant potential, in terms of them. From now we will 
call "A-representation" of the trilinear Higgs couplings their expression in terms of the Aj's. Moreover, 
when the Aj's are explicitly replaced by Higgs masses, v.e.v's and mixing angles, we will speak of "mass- 
representation". As far as the diphoton signal is concerned, the relevant couplings are those connecting the 
CP-even Higgs with the charged ones but also with the charged Goldstones. In the mass representation 
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they are given b; 



9hOH+H- 



m. 



vV^ { " V cos/3 
3 sin 2/3 



sin^2 COS/32 , flc , • flc ^ 

Jii + —r~7r'^i2 I + 2m^±(cos/3iii + sinfjbi2) 



sin/3 



m2 p'2 



+ 



5, 



i2 



cos/3 ' sin/3 




(3.4) 



9h^G+G- - 



v^f2 



m^o (cos/3 Sii + sin/35i2) + 2M^(5i| 



where m^o, m^o and nT-^± are the physical masses and the mixing elements 5*, 



^] 1 ^ij 



form the matrices 



diagonaUzing the effective mass matrices Eqs.( 1.20|1.22 ). Note also that in the non-linear gauge the 
h^G^G" couplings depend explicitly on the gauge through the parameters 6i. These parameters also 
appear within the ghost sector in the couplings h^c^c^, where c^,c^ are the charged ghost fields. The 
non-linear gauge parameter a also appears in the course of the calculation. It originates from couplings 
with physical fields like W^W^^, W^W^^^ and unphysical ones: G^W^j, c^c^7, 0^0^77 (see for 
example |51|). The latter quartic coupling emerges from the non-linear gauge condition only. All these 
couplings arise purely from the ghost and Goldstone part of the gauge sector and are not modified by 
the effective potential of the Higgs sector. 

For the numerical evaluation, as explained before, we obtain the Higgs masses and the mixing elements 
from NMSSMTools and the values are fed into SloopS through the SLHA interface. Here we use the 
second possibility offered by NMSSMTools, to compute the Higgs masses following |48j. There pole-mass 
corrections can be taken into account or not. In both cases the mixing matrices are computed in the 
effective potential approximation (i.e at = where p is the external momentum entering the Higgs 
self-energies) . 

As an illustration of the gauge invariance of the parameter reconstruction (Eqs.( A!l|A.2 )), we considered 
two benchmark points from [54] , named NMP2 and NMP5 after the conventions in [54j. Their respective 
Higgs sector parameters are recalled in Table [T| together with the soft SUSY breaking masses of the stop 
sector Mq3 , M^^ and of the gluino sector M3. All remaining soft masses and trilinear parameters that 



are not given in Table [T] are set at 1 TeV. The resulting Higgs spectrum is summarized in Table [2] and 



Parameter 


NMP2 


NMP5 


tan /3 


2 


3 


A 


0.6 


0.66 


K 


0.18 


0.12 


/Ueff. [GeV] 


200 


200 


[GeV] 


405 


650 


A, [GeV] 


-10 


-10 


Mqi [GeV] 


700 


600 


Mi'^ [GeV] 


700 


600 


M3 [GeV] 


600 


600 



Table 1: Benchmark points taken from [54] for the diphoton decay width. 



computed within NMSSMTools according to the procedure 48 with and without the pole-mass corrections 



In addition to the masses given in the SLHA output we also need the mixing elements Sij and P/- to 



obtain the couplings entering the diphoton decay width. In the effective potential approach used in 48 



^For the sake of clarity we reproduce the expression of the charged Higgs couphng. It can also be found in appendix [b} 
together with its general expression in the A-representation. 
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Mass 


NMP2 


NMP5 


[GeV] 


no pole 


pole 


no pole 


pole 




129.4 


126.5 


96.1 


95.6 




133.1 


132.4 


128.9 


126.5 




470.8 


464.5 


659.9 


655.8 




116.4 


115.7 


93.9 


93.2 


468.7 


462.8 


660.1 


656.5 




454.4 


454.5 


644.8 


644.9 



Table 2: Higgs spectrum of the benchmark points considered. For the NMP2 point the lightest CP-even 
Higgs is SM-like and for NMP5 it is the second-to-lightest one. 



that we henceforth denote as EPA, they are obtained by diagonalizing the radiatively corrected Higgs 
mass matrices in the DR-scheme anew. However the definition of the diagonalizing matrices is ambiguous 
since the self energies entering the radiatively corrected mass matrices depend on the external momentum, 
the rotations matrices are defined as those that diagonalize the mass matrices at 



In 
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0. Then, 

whether pole-mass corrections are taken into account or not leads to the same diagonalization matrices 5 
and P': This can be seen formally as a missing higher-order correction, but, in this fashion, the physical 
Higgs mass matrix would not correspond to the DR one in the EPA. This inconvenience is circumvented 
in the physical effective approach (PhA) as we force the radiatively corrected mass matrices to be the 
physical ones by imposing Eqs.( |A.l|A.2| ). 

To reproduce the EPA method using pole-masses within SloopS, we take the DR-masses for the recon- 
struction of the potential of Eq. (1.4) but the pole masses for the kinematics of the process. Therefore the 
Higgs mass appearing in the coupling of Eq.(3.4) is the DR-mass, which differs from the energy at which 
the decay is evaluated. This mismatch in the EPA with pole-masses will lead to a violation of gauge 
invariance within our generalized non-linear gauge, as we will show numerically. On the other hand, in 
the PhA, we reconstruct the potential of Eq.(1.4) (and consequently the Higgs mass matrices as well) 
directly from the pole-masses, which are also still used in the kinematics. This procedure will guarantee 
that gauge invariance is maintained because Eqs.( A.l|A.2 ) are fulfilled. At the numerical level we will 
vary the parameters a, 5i within SloopS to exemplify the gauge-invariance of the calculation in the PhA 
method. The results are presented in Table [3] They were cross-checked with the standard NMSSMTools 
version and with NMSSMTools*. The corresponding output of NMSSMTools, whether in the standard or 
modified version, does not check internally any consistency requirement, such as gauge-invariance, and 
simply uses an analytic, pre-computed expression for the effective /i^77 couplings. The excellent agree- 
ment among the results of NMSSMTools and the computation of SloopS for vanishing non-linear gauge 
parameters (i.e. in a linear gauge) is therefore a welcomed feature. The sources for possible discrepancies 
between SloopS and NMSSMTools lie in the treatment of the sfermion sector: in SloopS it is treated purely 
at tree-level whereas NMSSMTools includes several corrections to the spectrum and couplings 32 . For 



the two benchmark points investigated this difference is almost invisible because the sfermion sector is 
essentially decoupled. 

These illustrative examples show however evidently that the EPA calculation is not gauge-invariant. The 
origin of this breakdown can be traced back to the observation that Eqs.( A.1|A.2 ) are not satisfied. In 
more restricted gauges the gauge dependence would be seen at higher orders only. As stated above, 
setting a = — 1 removes the G^W^^ coupling and varying 5i^2 then gives gauge-invariant results in both 
procedures: they differ only by a finite and gauge-independent piece. This is due to the fact that the 
5i^2 gauge-dependent parts only appear proportionally to On the contrary, in the general 

case, d-dependent parts are proportional to the Higgs mass m^o, originating from the kinematics (i.e. 

i 

the center of mass energy \/s), and the g^oQ+Q- coupling, see Eq. (3.4). Recall that in the EPA the 
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NMP2 


r(/i? ^ 77) 


a = 5i = 


a = 5i = 10 


Sloops (EPA) 


1.138108952362.10^5 


4.490893854783. 10-5 


Sloops (PhA) 


1.125710969262.10-5 


1.125710969261.10-5 


NMSSMTools_3.2.0 


1.12699441.10-5 


NMSSMTools* 


1.12737737. 10-^ 


NMP5 


r(/iO ^ 77) 


a = 82 = 


a = ^2 = 10 


Sloops (EPA) 


1.053756232511.10-5 


3.628709516521.10-5 


SloopS(PhA) 


1.044860481657.10-5 


1.044860481613.10-5 


NMSSMTools_3.2.0 


1.04342526. 10-5 


NMSSMTools* 


1.04361857. 10-5 



Table 3: Gauge invariance test for the computation of T{h^2 ~^ 77) (i^i GeV) in the EPA and PhA 
procedures. Only the PhA approach passes the gauge invariance test within SloopS. There is no such 
test available with NMSSMTools, whether it is the modified version or not. 



mass appearing in the coupling g^^oQ+Q- is not equal to the pole-mass, which is used for the on-shell 
decay. The gauge-dependence of the EPA is precisely caused by this mismatch between the "kinematical" 



mass and the mass appearing in the coupling of Eq. (3.4) if pole- mass corrections are applied. In the 
opposite case, that of the PhA, the "kinematical" mass and the one appearing in the coupling are the 
same and the gauge dependent part oc d vanishes. We can render the EPA approach gauge-invariant if 
no pole-mass corrections are applied (this would then be a 'PhA with Di?-masses'). However, in this last 
case some precision is lost since, looking back to Table |2] there is a 2-3 GeV mass difference between pole 
and running masses. Within our reconstruction, the charged Higgs contribution is modified through the 
9h°H+H- coupling: we do not expect significant modifications with respect to previous calculations, in 
the MSSM-limit of the NMSSM, since it is known that the charged Higgs contribution to the diphoton 
decay width in the MSSM remains small (see for example (12]). In the NMSSM with large A (a form of 



the so-called 'A-SUSY' models 16 , which typically leads to a Landau pole below the GUT scale) and 



a relatively light charged Higgs mass, one could modify significantly 17^0 ^+^-, without requiring large 



doublet-singlet mixing. This was explicitly shown in 55 . 

Note however that a serious issue would arise with the gauge-dependent calculation of — t- 77 if one 
would choose to use it in order to derive some fundamental parameters at the Lagrangian level, which 
should preferably be determined from gauge- independent observables. As a final remark concerning this 
section, beyond maintaining gauge invariance, the PhA is clearly advantageous as it enables us to use 
the pole-masses easily in the calculation of the decay width, in a consistent way, without resorting to the 
technical task of computing the "pole-corrected" mixing elements Sij , P/^ . 

3.4 Comparison with micrOMEGAs 

A comparable approach, based on an effective potential approach, had been carried out in micrOMEGAs 
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34j, a code computing the DM relic density ^^h? in (e.g.) the NMSSM. Since light Higgs states can 
be present, annihilation channels into h^h^, h^a^, a^a? can contribute significantly to and such 

channels are affected by radiative corrections in the Higgs sector. The effective scalar potential was 
implemented in this code as (see |34|, where a slightly different version was proposed), 

Had = Xf\Huf + Xl'm^ + \i'\H^\^m^ + Xf\Hu.Hdf + ^[{H^.Hdf + h.c\ 



+ Xl\S\^\H^\^ + X'2\S\^\Hd\^ + 



^S\Hu-Hd) + X;S^ + h.c 



+ yl^r (3.5) 
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Parameter 


Value 


Parameter 


Value 


Ml [GeV] 


84.49 




2 


Ms [GeV] 


359 


A 


0.63 


M3 [GeV] 


1200 


K 


0.05 


Af [GeV] 


-1500 


Ax [GeV] 


694 


Mj [GeV] 


200 


[GeV] 





Mq [GeV] 


600 


^eff. [GeV] 


300 



Table 4: SUSY point for the comparison between micrOMEGAs , SloopS and NMSSMTools. and Mg 
are common sleptons and squarks soft masses. 



This potential is to be understood as a radiative potential, which means that all the parameters are 
loop-induced. To make a connection with our conventions we have 



AAi 



Af , AA2 



^M 



AA, 



Af , AA4 = Af , AA^ = XI , AA^ = A| , Ak' 



Af 



(3.6) 



where the AAj are the loop-induced part of the Aj parameters appearing in Eq.(1.4) when we split them 
as Aj = A^-|-AAj. The remaining parameters A^, Af , and A5 have no equivalent in our restricted potential 
of Eq. (1.4), but correspond, in our conventions for the general potential (Eq, 1.2), to 



^5 



Ar 



Ap 



4/4 



(3.7) 



Conversely, our parameters A^d, As and A 
obviously, does not rely on the 



M 



15 ana A'p receive no correction in the micrOMEGAs approach, which, 
^3 symmetry. As a consequence, while an inversion procedure is also 
possible with the potential of Eq. (3.5), the radiative corrections to the masses will be distributed in a 



different way among the Aj's, leading to differences at the level of the Higgs self-couplings. 

We remind here, that if one aims at improving on the tree-level couplings, as is obviously the purpose 



of a radiative potential, i.e. of Eq.(3.5), it becomes crucial to identify the Aj's that are subject to large 



quantum corrections: that was our discussion in subsection 2.3 An arbitrary truncation of the potential, 
albeit allowing for an inversion in terms of the Higgs masses provided it is sufficiently simple, is not 
a receivable option because the accuracy contained within the couplings brings no improvement with 
respect to the tree-level evaluation. Our study of the large logarithms within the Coleman- Weinberg 
approach tends to convince us that our choice of a Za-invariant potential should be preferred, while the 
choice in Eq.(3.5) seems arbitrary. Possible reasons for this choice within micrOMEGAs could lie on the 
facts that the loop-corrections to A5 are sizeable in the MSSM, and one could have expected the same 
behavior in the NMSSM, while the parameters A* and A5 (or rather k| and Af in our conventions) 
appear in the trilinear couplin gs h^h^h^ and h^a^aj. with a factor s (questionably an enhancement factor 
in the MSSM limit): see Eqs.( B.5|B.6 ). However other parameters in the general potential of Eq.(1.2) 
share this latter property and are still arbitrarily absent. We will see in a numerical example that ensuing 
deviations between our implementation and that in micrOMEGAs could be significant. 
We considered a point in the NMSSM parameter space where the DM relic density Q^h"^ is in the correct 



experimental range (at the 2a level: 0.1 < ~ 0-124 M), when computed with micrOMEGAs_2 .4.1 



and the Higgs radiative potential of Eq. ( 3.5 ) . This specific point passes warnings from NMSSMTools_3 .2.0 



as well, and features a SM-like CP-even Higgs mass around 125 GeV. The NMSSMTools input for this point 
is given in Table |4] The main channel contributing to i^-^/i^ is XiXi ~^ a^a^ (at 72%) and the rest of the 
contributions involve fermions in the final state, dominated by the bb final state. The process XiXi ~^ aifl? 
is dominated by the s-channel exchange of the SM-like Higgs close to its mass shell, as can be seen on 
Table 5| where the Higgs spectrum, the lightest neutralino mass m^o and the resulting relic density are 
provided. 

The decay modes /12 3 ~^ ^a^ are kinematically open. Let us compare these decay widths within the 
three codes: SloopS, with our effective implementation (Eq.( A.1|A.2 )), micrOMEGAs using the radiative 
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Spectrum 


m^o [GeV] 


63.2 


m^o [GeV] 


110.9 


m^o [GeV] 


126.4 


m^o [GeV] 


727.8 


m,o [GeV] 


59.7 


m,o [GeV] 


732.4 


m^± [GeV] 


721.7 


Relic density 




0.103 



Table 5: Resulting spectrum and Q^h'^ from the data point presented in Table [4] 



potential in Eq. (3.5 ) and NMSSMTools_3 .2.0, where only the leading logarithms in top/bottom corrections 
are taken into account. The output of NMSSMTools* is also considered. The results are displayed in Table 
rol We observe significant discrepancies between SloopS/NMSSMTools*_3.2.0 (which are in remarkable 



Decay [GeV] 


m ^ «?«?) 




SloopS 
micrOMEGAs_2.4. 1 
NMSSMTools_3.2.0 
NMSSMTools*_3.2.0 


3.566 10"^ 
2.96010-2 
2.73010-2 
3.56610-2 


1.90010-4 
4.66510-^ 
1.23310-4 
1.90010-4 



Table 6: Comparison of the decay widths 3 ~^ ^a?- 

agreement), on one side, inicrOMEGAS_2.4. l/NMSSMTools_3.2.0 (which also show some disagreement 
between them), on the other. Our calculation for the main channel h2 — )• 2a^ (that we denote henceforth 
as T^{h2 alal)) is about a factor 1.2 larger than the micrOMEGAs result (labeled as r^(/i^ «iai))- 
Giving the modified prediction of O.^h'^ within our procedure is beyond the scope of this work, but we 
can nevertheless make a rough estimation of this quantity. As the process XiXi ~^ a^a^ is dominated by 
the /i^ resonance, and only the h2aiai coupling is modified, we can reasonably approximate, 

a^(x?X? ^ aK) r^(/.° ^ gggg) 

tT^(x?X? ^ a?a?) " rA^(/i^ ^ a?a?) " 
Denoting cr*^ as the contribution of the rest of the processes to the cross- sections involved in Q^^h"^, the 
relic density computed within micrOMEGAs, we can write the sum of all contributions 0"(|^,-(xiXi ~^ ^) ^s, 
X standing for any relevant final state, 

-t'^t(x?X? ^X)=Y: <r'\x\x\ ^X) = a^^(x?X? ^ a?a?) + a^" (3.9) 

X 

The micrOMEGAs calculation gives 0"*^(xiXi ~^ Oi'^i)/<^tot(XiXi ~^ X) = 72%, as we already men- 
tioned, and the ratio of relic densities in both approaches is approximately determined by Q.^h? /Q.^h? ~ 
^tot/'^tot) with Q.^h'^ the relic abundance in our calculation, and where 

Moreover we have = a^'^ , since the remaining relevant contributions are annihilations into light 
fermions and hence unaffected by corrections in the Higgs sector. Thus we obtain the following estimate, 

n^h'^ ~ 0.090 (3.11) 
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A reduction of the relic density with respect to the micrOMEGAs calculation was to be anticipated since 
in our computation the annihilation into light pseudoscalars is enhanced, thus depleting the abundance 
of relic neutralinos more efficiently. In turn, and contrarily to the prediction of micrOMEGAs, this point 
would actually lie outside the cosmologically interesting region if one relies on our estimate. Of course the 
derived value of ^x^"^ depends crucially on the precision of the evaluation of rUf^o (and m^o) since the 
annihilation XiXi ~^ a? a? occurs at the /ig resonance. These considerations are of particular significance 
when one considers that the PLANCK satellite |57| should improve the experimental determination of 
cosmological parameters |58| soon. For a discussion concerning the accuracies required from colliders to 



match the precision of the relic density measurement, see for example 59 



Conclusions 

This study of the Higgs potential with two Higgs doublets and one gauge singlet has put forward several 
points of interests that we would like to summarize briefly here. 

The most general effective renormalizable Lagrangian of the 2-doublet+l-singlet setup, contains 28 (plus 
one superfluous) parameters, far beyond the 11 ones of the 2HDM, even after complex phases have been 
discarded. Therefore, if future experimental measurements should point towards such a rich Higgs sector, 
a full reconstruction of its potential through experimental data in the Higgs sector could succeed only 
after an exhaustive measurement of the Higgs self-couplings, beyond that of the masses and mixing 
angles: if the purpose for such a reconstruction is sound from the point of view of model identification 
and precision tests, it is also probably condemned to a very long delay, as far as the experimental phase 
is concerned. This situation is eventually that of most models, albeit constrained, once considered at the 
radiative level, for symmetries are spontaneously broken by the Higgs v.e.v.'s and loop corrections end up 
contributing to all possible terms in the potential. We emphasize, however, that a precise determination, 
in a general parametrical form, of the potential at future (linear) colliders, shall help discriminate among 
such models and constrain their parameters: in turn, the predictions of specific models for the parameters 
of the Higgs potential should be known at the radiative level so as to allow for comparison/precision tests. 
Requirements for additional symmetries, beyond the EW-invariance, or matching conditions originating 
from more-elaborate models may constrain the effective potential at the classical order. Provided its form 
is simple enough, an identification at leading order of the parameters of the underlying model is achievable 
from the Higgs spectrum solely. Then, assuming the remaining sectors of the model are sufficiently 
documented as well, a full determination of the effective potential within the more-fundamental model is 
essentially a matter of perturbative calculation. We lent particular attention to the Zs-invariant and PQ- 
conserving potentials, which could both be embedded within a SUSY extension of the SM, respectively, 
the NMSSM or the UMSSM, the PQ '-conserving potential (R-symmetric limit of the NMSSM) or the 
potential driven by an underlying nMSSM. A reconstruction of the classical parameters was explicitly 
carried out, at leading order, for those models. 

Further achievements seemed within reach in models ensuring a residual symmetry at the EW scale. Our 
test-model here was the NMSSM, and the study of the large logarithms within the Coleman- Weinberg 
approach confirmed that the leading-logarithmic effects would not spoil the Za-symmetry, extending the 
validity of our parameter reconstruction in terms of the Higgs spectrum to this order. By contrast, in the 
nMSSM, where no residual symmetry is present at low-energy, logarithms do not observe the classical 
form, spoiling a reconstruction beyond LO. 

We finally considered a few phenomenological consequences of this parameter reconstruction at the 
leading-logarithmic order in the NMSSM. We based our discussion on the Higgs spectrum computed 
in the public code NMSSMTools and implemented the reconstruction both within NMSSMTools, directly, 
and within SloopS. The latter allowed us to visit the diphoton decay of the SM-like CP-even scalar again 
and clarified the conditions for a gauge-independent implementation. Comparison with the previous 
implementation of an effective Higgs potential within micrOMEGAs was also carried out: different choices 
in the radiative potential result in different Higgs-to-Higgs couplings at the order of leading logarithms, 
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as the radiative effects encoded witliin tlie masses are distributed differently among the parameters of 
Eq.(1.2); while the form in Eq.(3.5) is seemingly arbitrary, our choice (Eq.(1.4)) is justified by the anal- 
ysis of the logarithms appearing in the Coleman- Weinberg approach, and should thence prove a priori 
more reliable. As far as the phenomenology of the NMSSM is concerned, we found fine effects in collider- 
constraints or the calculation of the DM relic density, appearing essentially for points of the parameter 



a^a? or 



space which rely heavily on Higgs-to-Higgs couplings, such as those entering the processes 
XiXi ~^ ^i^i, mediated by a CP-even Higgs in the s-channel. 

Finally, let us mention that, although the state discovered at LHC is in a favourable mass-range for singlet- 
extensions of the MSSM, a long stage of experimental measurements and identifications of additional 
Higgs states lies ahead of us, should the 2-doublet-|-l-singlet setup be realized at all in Nature. 
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A Parameter reconstruction for simple classical potentials 



We provide here the results of the inversion procedure described in section 2.3 for a few classical potentials. 
Note that, for completeness, one should also replace the parameters within Eq. (1.15) to fully determine 
the potential. 

Za-invariant potential V^^: 
The quartic doublet couplings are entirely determined by the Higgs mass-matrices: 



Ai 

A2 

A3 
A4 



1 

2v^ 



1 



1 
2v^ 



1 



rrp- P'2 



sin^ 13 tan2 /3 

2 D/2 ^2 



2m2 ± + 



(A.1) 



One degree of freedom remains, which can be chosen conveniently as the singlet v.e.v. s. One then obtains 
for the remaining parameters: 



Ap - 



1 
4s2 



s av ' ^-1- 

sin 2/3 2 p/2 _ 12 pi pi 

,,2 „;„2 



^ sin 2/3 2 p/2 p/i! "sinzp y p/ p/ 

2i2 "^aO-^il ~ ^a°^i2 2^^oP^il^i2 



,2 Z3'2 _ t^sin 2)9 ^2 p/ p/ 
i2 

,2 c2 I lrrj2 p'2 sin^ 2/3 2 p/2 



(A.2) 



Ap 



i 

2sv sin /3 

i 

2sv COS /3 



+ 



+ 



1 

3s tan /3 



sin 2/3 2 p/2 _ 2 pi pi 



tan / 



sin2/3^2 p/2 1 



^2 pl2 1^ -Z pi pi 



In replacement of s, one may use any combination of these latter equations to define a new parameter. 
For instance, 



-ia,b) 



aA^ + bXj, 



44> 



V sin 2 l3mloP^^ 


[l 2e(,,,)(j^^+6tan/3)] 




9^2 pi pi 

a" *1 j2 

i 




+ 3e(,,,)m2„5,3 (« + ^of/3 ) 



(A.3) 



£(a,b) coincides with k/X in the NMSSM at tree-level and may be regarded as a measurement of the 
breakdown of the Peccei-Quinn symmetry. Alternatively, 



Pud 



Ap 



2v sin 2/3 mloPlf [l - (1 + Pud) tan^ /3] 



^loPkPL [1 - (1 + P^d) tan2 /3] - 3m2o5i3 



cos B 



'1 + /5«d)tan2/3 



sin /3 



(A.4) 



Pud vanishes at tree-level in the NMSSM and may represent another possibility. 
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Peccei-Quinn-invariant potential Vpg: 

The system is fully determined by the mass matrices: 



1 

2v^ 



1 

2v^ 



Ai = 

A2 = 

A3 = 
A4 = 



s = f sin 2/3^ 

t)2 sin2 2/3 



COS 2 /3 


2 p/2 


tan 


i 


"2 




sin^ /3 


tan^ /3 




2m^± + 


2m2QSii. 




sin 2/3 
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1 



2 p/2 

™2 p/ p/ 

"'0-^21-^22 

°2 



ra 



2 c-2 



2 p/2 

"^ar 22 



Ap 



Ap 





1 


p' 


m^QSi2Si3 

i 


2t>2 


sin 


PP2I 


COS /3 




1 


P' 




~ 2v^ 


cos^ 


PP2I 


sin /3 



+ "^^0-f21-^22 



"^„0-f21-f22 



(A.5) 



Peccei-Quinn'-invariant potential Vpg,: 
The system is fully determined by the mass matrices: 



Ai 

A2 

A3 
A4 



Ap 



1 

2t)2 



1 

2v^ 



1 

2i^ 



t 

COS 2 /3 



sin2 /3 



m^„P^^itan^/3 



"2 
tan^ /3 



2m 



+ 



2m.2gSiiSi2 
sin 2/3 



Ap 



Ap 



2 p/2 

- °2 
2»;2 sin 2/3 

sin 2/3 §1 

^22 

4t)2 sin2 2/3 V 7^ 



41)2 sin2 ^ p/^ 



4?)2 cos2 /3 P^^ 



"^00-^21-^22 



"^aO -^21-^22 



m2 



m2g5i25i3 

i 

cos /3 



i 

sin /3 



(A.6) 



nMSSM-inspired potential V|^: 
Although only twelve parameters are to be determined within the potential, application of the constraints 



of Eq. (|2.5|) leave one degree of freedom, due to the degenerescence of the CP-even and CP-odd singlet in 
2 02 



this model: m^nS}^ = m?^QP'^. We again choose s to be this degree of freedom. The remaining parameters 
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read: 



^1 - 2^^ 



A2 = 



2v^ 



^3 - 2^ 



cos^ p a" ^ 



,2 P'2' 



sm2/3 tan2/3 

2m2„± + P'2 

i/± sin2p a" «i 



i 



\a — — 5" 

"^?2 = -^oPA^^2f-rn^oP/?#^ 



(A.7) 



m 



I pi pi 



2vs 



\d _ 1 
■^P — 2vs 



V 

m?,,Si2Si-i m? 



pi pi 



sin /3 



tan /3 



B Trilinear Higgs-to-Higgs couplings 



In this appendix we give the physical trihnear Higgs-to-Higgs coupHngs h^H^H , h^oP-a^ and h^h?-h^ in 



i j' k 



the A-representation obtained from the general potential Eq.(1.2) and in the mass representation from 
the restricted Z3 potential Eq.(1.4) only (as in the general potential the results are cumbersome). In the 
following the matrix Pij is defined as the 3x3 diagonalization matrix which rotates the gauge eigenstates 
(0^,0^,0^) directly to the physical basis {a^,a2,G^) such that, 

■P[,s^ P[,cp P[, 

P^l^P -^21^/3 P22 
C/3 -SjS _ 







0.2 









Pu 


P12 


Pl3 


P21 


P22 


P23 


P31 


P32 


P33 




(B.l) 



with P' defined from Eq. ( 1.20 ). To cast the couplings in a more compact form we also define the following 

(B.2) 



mixing elements combinations, 

Sia [PjbPkc "I" PjcPkb] 



(B.3) 



B.l Trilinear couplings in the A-representation 
CP-even Higgs to charged Higgses h^H^ coupling 

9hOH+H- = ^ \ ^T^r h V2A3t;[c^S'ii s^S'i2] 



V2 



V2 



[s2l3Sil 



V2 

(1 - 3c2^)Si2] + [(1 + 3c2p)Sa 



V2 



(A4 + X5)VS2I3 

- 8213812] 



[sibSh + Ci3Si2] 



+ 
+ 



S2P 

V2 



Aud + + 4(A|f + A^-' + \m)s) 



:m 



Si3 



Ads + (A^ + 2A^)s )sj+(A^s + (A^ + 2A^)s ) c} 



S, 



i3 



(B.4) 
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CP-even Higgs to 2 CP-odd Higgs /i°a°a° 



9h0oP.al 

I J k 



+ 



+ 



+ 



-^l^Cff ^^1^1,1 ^ A2t^,jj^. 2,2,2 ^ (A3 + A4)t; 

^2 ''^'^ V2 ''^''^ ^/2 



^t^^/ttAn 1,1,1 
V2 ^ ^^'•'■''^ 



^/2 



Xevc0 ^^ 2,2,2 _ ^^(Asc^^^Aes^ ^^ 1^2,2 _ -t;(A5S/3 - A7C/3) ^^ 2,1,1 



V2 



V2 



- V21 



An1,2,1 



^ud — 2Ap S ,^.1,2,3 Aid — 2ApS^ 2,1,3 I A«H-2ApS . 3 i_2 



V2 



V2 



^Mcf — 2A;^s ^^^^1^2,3 — 2Ap^s 2,1,3 , + 2(Ap^ + AAf)s ^ 3^1 2 



%/2 



V2 



-(n") 



+ 



V2 



-(n^ 



(B.5) 



Arf, + 2A^s . 3,1,1 + 2A|,s . 3,2,2 3^5 - ^5 - 3(2^2 _ 3^2 )s 



(n 



V2 ^" '^'^'^^ ^ ^/2 '^'^■''^ 



2(A|f-A|f). .^(n-);';^,^ + c,(n-) 



A\3,l,3 



TA^3,2,3 



+ 



V2 



3V2 

■^pS /i-rA\3,l,l , ^Pf m^^3,2,2 



3,3,3 

k 



V2 



Triple CP-even Higgs coupling h'^h^h^ 



9h°h°h0 

I J k 



+ 



>i,j,k + ^ 



li,j,k ' ^ 



^l!^m5^i,ia , ^2vs^. s. 2,2,2 , (A3 + X4)v 

XqV_ 

V2 L' 



/Tr5\2,2,2 , g /tt5'\1,2,2' 



'ijk 



X-jV 



c.(n^)^^:f + .,(n^)JS 

2,1, 



^Mn^)S + 3c,(n-)^'^'^ 



(A^c/3 - Af 5 1,3,3 (A;^S/3 - Af c/3)t; 5 2,3,3 , A' 



V2 



(n 



'i,j,k 



+ 



V2 



n^) 



(n 



5n3,1,1 



Xls 



(n^) 



5^3,2,2 



{2Xj,Cp-{XM+~Xj^)sp)v ^^Sp3,3 ^ {'^~X^pSf}-{XM + ~X^p')cfl)v ^^S^2,3,3 



^/2 



' V2 



S\3,l,2 



+ 



Sn3,3,3 



^„d + 2(AM + A|f)s 5 3,1,2 , ^5 + 3(4 + 4)s 



(n 



'i,j,k 



+ 



3\/2 



■(n 



5^3,3,3 



B.2 Trilinear couplings in the mass-representation for the Z3-conserving potential 

To obtain the mass representation we trade the Z3-conserving Aj's of the couphngs in the A-representation 
(see previous subsection) against the masses, mixing angles and v.e.v.'s (Eq.( A.1|A.2 )) and set the re- 
maining ones to zero. 



CP-even Higgs to charged Higgses h^H^H coupling 

1 



9h°H+H- 



j\/2 l"*^'"' V cos/3 
2ml,P'^^P'^2Si3 
3 sin 2^ 



sin/32 cos/32 \ „ 2 / «e , • fle 



sin/3 



Sil ^ Si2 



4 V 



"^aO-P/l^ \ a ' ■ Q o 
J ■' \ cos/3 sm/3 6 s 



Si3 



(B.7) 
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CP-even Higgs to 2 CP-odd Higgs h^a^aj 



2 

hO 



^ l-i 2^ 



1=1 ^ 



c2 



+ Sis 



+ 



+ 



+ 



02 

tAn3,1,1 



+ + + ^^^^ 



5n 



1.^^ . )i,j,k \ , c._ / 



PZ1PZ2 



/j^A\l,2,2 



+ 

VCp s 



2vc 



■/3 



+ 



+ 



x2,3 



+ 



(n^)2,l,l 



VCjj 



2vs 



•0 



1 / (n-^)-g , (^tii (ji^tii 



+ PllPl3 



1 



(n 



3s 

+ P12P13 



2v 



A\?>,1,?> 



3s/= 



l/-rrA\3,2,3 



fn 



Axl,3,3 



+ 



2s 



1 

2^ 



C/3 



+ 



C0 



■/3 



3s 



1 



3s 
2s 



(n 



Aa3,2,3 , , /T-rA\3,l,3 /t-tAn2,3,3 ^P^^ 



+ - 



A^3,3,3 



2s 



(B.8) 



Triple CP-even Higgs coupling h^hP-h\ 



+ 



m. 



^2a/2 1 wc/j i;s/3 



1=1 



+SiiSi^ 



+ 



VC0 
+ 



+ 



(nS)l,2,2' 



'i,j,k 
VSfi 



US/3 



+ 



n3,3,3 
'i,j,k 



p2 
Ml 



UC/3 



+ 



rSN3,l,l 



lm5'\3,l,2 



+ Pj 



12 



2 //TT5^1.3,3 /TTS^3,3,3 



2s0 



6s 



(B.9) 



r5x3,2,2 



vs 



nS'\2,3,3 



-P<l-P<3 

3 

P12P13 



4 

3s '''^'^ ■ "'^'^ ' s \ 2cp 

^ij.fc ^ _ V-^^ h,j,k ^ VCp .-^5x3,3,3 

UC/3 VS^ S 6s2 

("^)Sfc ^ ("^)Sfc _ ^ !!f^m^)3,3,3 

6s2 '^d^k 



/'tt5\3,3,3 

6s 



i;s/3 



+ 



^'?3(n^ 



'i,j,k 



9s 
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C Coleman- Weinberg analysis of the Higgs potential in the NMSSM 
and the nMSSM 



The two models under consideration essentially differ, at tree-level, by their Higgs sectors. Additionally, 
one should require the limit — in the nMSSM neutralino sector, with respect to that of the NMSSM. 

SM-fermion contributions: 
In the base of Dirac-fermions (u,d,iye,e), the squared mass-matrix of SM-fermions in terms of neutral 
Higgs fields reads (we omit color and generation indices) : 

r0|2 



leading to the potential: 



u \ u\ 
















































Y2 

^ e 




2 


In 




^) 


^\ 


A2 


) 



(C.l) 



A2 



+ 



(C.2) 



where we have kept only the leading, S'?7(2)i-invariant, logarithmic terms. We deduce the corresponding 
contributions to the Higgs potential: 



a/ ^ V y4 1 / 



(C.3) 



SM-Gauge-boson contributions: 

In the base of real vector fields (7°, ^ W2, Z°), the squared mass-matrix of SM-Gauge-bosons in terms 
of neutral Higgs fields reads: 



1 



■ 
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^2 Q 

g'^ + g''^. 



(C.4) 



leading to the potential (note that in the S'C/(2)/,-conscrving limit, these fields are massless): 



5V^^{HlHlS) 



647r2 



2 



+ 



(|^«| +\^d\ 



In 



In 



2A2 



'{g^+g'^)[\Hl\' + \Hl\')\ 3' 



2567r2 

providing us with the couplings: 



2/ In 



A2 



(ff^+5'^)^ln 



Ml 
A2 



2A2 

(l^^l 



(C.5) 



T287? 



T287? 



A? + X9 




+ ig^+g'^fln 
+ {g^+g'^ fin 



'Ml 



Ml' 



(C.6) 



1287r2 



+ {g^+g''^f\n 



Ml 
A2 
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Sfermion contributions: 

The Sfermion squared mass-matrix, in the base (Fl, F^*) for a flavour /, is given by: 



MUS,H"^,Hli) 



m% +\YfH<}\^ + 



m'^-m'') Yf{AfH'}-XS*Hj*) 



Defining = (A^|)n + (A^|)22 and i?| = (A^^ )n 



and the Higgs potential: 



(A^|)22 



(C.7) 

^2) we obtain the eigenvalues 



efF v"«>"rt>-/ -^287r2 




+ 2rp -i^^ln 




(C.8) 



One derives the couplings: 



fid,e)Yf 



2r/ 



9"y 



Al+Ai + A^Bf,, 



4 ■^L 



1671-2 



/ |"/(<i,e) 



2y/ 



g'\yl+yi,)-2g''ii 



In 



A2 



In 



A2 



Y2A2 a (yj-yD-^a H 



In 



A2 



In 



-^2 



\Yf + 



4 



,2 r/ 



In 



A2 



A2 



In 



A2 



-YfA} 



+6 



fid.e) 



2YJ 

9'-yl-2g-li 



g"{yl+y'n)-2g'Ii 
4 

2 

In 



4y 



g"(yi+y;,)-2g^/( 



A2 



In 



9l01zZiilz^C£iA^2i? 

4 ^"''L,R 



r)F 



(A3 + A4)^ = T^E/ 



i(r/ 



■-LM 

g'^yi-2g 



A2 



A™2F 
4 ^^L.R. 



In 



A2 



g'^y[-2g^lf 



Yf + 



In 



A2 



+y/AV[(F/ + ^:^M±f)z^ 



^(l+Afe))cf,H-(24 

i(y/ 



+ 



(l-5^)Ami5j]7?r,^ 
In 



g'2y/-2g2j/ 



A2 



+y/A2s2 



Af 



167r2 ^f^J 



Y^AfXs 



(A 



u \F 



+Sft^a,e)Y^AfXs 

6fi^d,e)YfAfXs 
+SfuY^AfXs 

1 V v2\2rF 
1671-2 /^f=d,e f '-'L,R 



g'^(y[+y{,)-2g^lf ^F 



16^2 E/ { 



■'L,R 

g'^(yi+yi)-2g^ii 



cf,;«-(i7AV + 



2^2„2 _L g"(yi-y^)-2g^-f| 



g'^{y[+yi)-2g^li 

8 '-'i.-R 



^L,R 



g'^(y^-yD-2g2j/ 



2A2s2 



g'-(y^-yD-2g2j/ 



. {^py — li^ E/=« Yf^^^LR 



(y/ + --'^^"^^^^---'^^^ (i + ^)) cf,, - (24 + ^'^^M.^l^l^ - i)A<,) p£ J] } 



(C.9) 
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where we have used the notations: 



Af. = ml 



In 
1 



A2 



L.R 



'^L,R — 



2 2 



In 



Ar = 



In 
1 



— 



In 

2 

m 
A2" 



A2 



In 



me, 



In 



A 2 F _ 2 2 



6m4 



(C.IO) 



Chargino contributions: 

The chargino squared mass-matrix, in a base of Dirac (winos,higgsinos), is given by: 



Ml^{S,HlH^,) 



(C.ll) 



Defining T^± = 
eigenvalues tz^^-i 



(>i'±)ii + (A^^±)22 and Rl 



One can derive the coupHngs: 




(A^L)22 



+ 4|A^J±I!2, 



we obtain the 



and the Higgs potential: 



+ RU) 



In 



'rp2 



4A4 



+ 2T^± • i?^± In 



(C.12) 



A 



ud 



(A 



Ml 



, (A«,<i) 




(Al,2)^ = 

(A3 + A4) = -8„2(m|-A2s2)3 



-AS6ln(4#)} 



16Tr2(M|-A2s2)3 



167r2(M|-A2s2)3 



2M| 



87r2(M|-A2s2 

+3M|AS4 



167r2(M2-A2s2)3 



2M| 



ln(f)- 
ln(4i^) 



ln(|l)+21n(4#) 
ln(4i^)+l]} 



+ 4M|A2s2 - AS^ 



(C.13) 
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Neutralino contributions: 

The (hermitian) neutralino squared mass-matrix is determined by its entries in the base of Weyl spinors 



{Ml. 
{Ml. 

iK" 

{Ml. 

{Mlo 

{Ml. 

{Ml. 
{Ml. 
{Ml. 

{Ml. 

{Ml. 
{Ml. 

{Ml. 



11 



22 
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<? ffO TirO 



^7 H^) 



5', -H^Si-f^d) 



S,Hli,H^,)=M! + '-{\H°^ 



0,2 , |„0|2 



99 



{w+m 



S,H',,H',) = X'{\S\' + \H^,n + 



0\2\ , 9^ + 9'^ \ JjO\2 



S, Hi Hi) = X\\S\' + \Hl') + ^--^\HX 



X'^S*H°-2XkSH^' 



S,HlH'i) = X'S*H'i-2XKSH]i 



0* 



^{M,Hl*+XSHl) 

-l={M,Hl*+XSHl) 
-j^{M,H^^*+XSHl) 
-j={M,Hr + XSHl) 



(C.14) 



One can expand its eigenvalues in terms of doublet fields: 



mj{S,HlH^) = m]^'\S) + mi^'\S,HlH^)+my{S,HlH',) + 0{{Hl,f) 



2(1) 



2(2), 



2 {2d) 



2 (Id) 
i 

2{2u) 



{S)\Hl>\'^ + [m^'^'-'\S)H^^Hlj + h.c. 

2 (2nd) , 



(C.15) 



The associated potential is then given by: 



SV^/ {Hi HlS) = ~^Y. ^t{S, Hi Hi) 



=-—y 



'{S) 



In 



mKS^HlHl) 
A2 



In 
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'2m^^'\s)m^^'\S,HlHl) 



A2 ; 2 



+ 2m^^'\s)m^''\s,HlHl) 



In 



m^''\S) 
A2 



In 



A2 



mt^'\S,HlHl)ln 



A2 



(C.16) 



Mlo{S,Hl^ = 0) is already diagonal with the eigenstates {Mf, l^^i)}, {M|, \E2)}, {X^\S\^, , {E^)} 
and {4«;^|5'p, l^s)} stands for the canonical base of C^), from which one obtains easily the 
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pure-singlet parameters: 



As 

2\X 



= -4^(A4 + 8^V 



-l6^{A^[ln(^) 



(C.17) 



+ 



The 0{H'^) masses also come without much effort: 
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2(1) 



2(1) 



- A2|5|2 



M| -A2|S|2 
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m 



2(1) 



+ m; 



2(1) 
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+ 



2A2 



(A2-4k2)|5|2 



2(1) 



8A2 



(4/^2 - A2)|5F + 1^^°!') - «ARe(5*2ij0if0)] 

from which we can derive the couplings (we focus on the logarithmic terms) : 
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- A^s^ [3M2 + A2s2] In | 
[Ml +3A2s2] In (^) 



A2 



A^ 
A2 



g^AMa 



Ml 



5A2 




A2s2 [2M4 



^ 167r2(Af2_A2s2)3 

[^uaf = - le^^^fSw + 3A^.^] m (^) - 



3M2a2s2 - A4.s4] In 
A2s2 [2M| + 3M|A2s2 _ A4,s4 



(a-.)' 



g'2A2 

■ 167r2(M2-A2s2)3 

9^A2 



167r2(M2-A2s2)3 



+ 



8ii-2(4k2_a2) 

/2 J, t4 \4 2 

g JV/^ A s 



167r2(M2-A2s2)3 



g'2Mj'A'^^^ 

167r2(M2-A2s2)3 



{m4[m2 + A2s2] In 
{m| [Ml + A2s2] In 
- 16k^ In 



kA-' 
47r2(4K2_A2) 

g'^M^'A-^a 



87r2(M2-A2s2)3 




A2 

M. 
IP' 

4k 



^)} 



^)} 
(^)} 

[Ml + 3AV] In(^) 
A2s2 [4Mf - 3M2a2s2 + A^s^] In (4§i) } 
A2s2 [4M| - 3M|A2s2 + A^s^] in (^f!) | 



g^MlA'^s^ 
167r2(M|-A2s2)3 



(Am)'^ 
^A^'l'' 

l,\^^My ~ 167r2(M2-A2s2) 

The limit -)• for the nMSSM is straightforward. 



ln(#) 



_ g2M|A*,s^ , 
167r2(M2-A2s2)3 
g2M|A"s3 

167r2(M2-A2s2)3 

4k2 In 

g^M|A^ 



_M£_ 

A2s2 



\A2s2 j 



87r2(M2-A2s2)3 

g^M|A^^ 
^ 167r2 (Ml -A2s2) 
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Charged-Higgs contributions — NMSSM: 

In the base {H~,H^), the hermitian squared mass-matrix of Charged-Higgs bosons in terms of neutral 
Higgs fields reads (we use the general notation of a Z3-conserving potential; those parameters should 
be replaced, in practice, by their tree-level value; we also define Ap, replacing Ap*^, which coincide at 
tree- level; same thing for A., replacing Ai,2): 



tan / 



{M%i)22{S, H%) = Ml^ + (Ad + A^s)stan^ + AJF^I^ + M\H'i\^ + Ap(|5|2 
(A^^±)i2(5,<,) = A^aS + X^S*^ - Xa{KhI)* 



(C.20) 



We have introduced M'^^^ to replace constant terms generated by the electroweak v.e.v. and regularizing 
the (otherwise- vanishing) Goldstone mass (which does not correspond to a Goldstone boson since SU (2) l 

is conserved in our approach). In practice, this M'^yy should be chosen as M'^^i typically, since it 
replaces the longitudinal component of I^-bosons. Now, defining T = TtA4^±{S, ^) and = 
- 4detA^|^± {S, H!^^^), we obtain the two eigenvalues m\i^{S, H^^^) = ^[T - / + R], as well as the 
potential: 



SV^f\HlH'„S) = 



1287r2 



(T^ + R^ 



In 



R^ 



4A4 



- 3 



+ 2r- i?ln 



T + R 
T-R 



(C.21) 



Then we focus on the logarithms (the notation (•) means that Higgs fields are replaced by their v.e.v. 's): 




In 



~ In 



1 + 



sin 2/3 
-I- . . . = - In 



+ . . . = In 



A4 



+ ... 



4- ... 



(C.22) 



Expanding their coefficients, we obtain the leading charged-Higgs contributions to the Higgs-potcntial 



parameters. Note that the coefficients multiplying In 



are in general very complicated. Here, for 



simplicity, we give only the leading term in sin 2/3 — > (tan/? — > oo) 



327r2(A„d+Ajfs) 



(5A„d + 12A^s)ln 



{^ly^ ^ 32^ { e^^sInS^ (^.^ COS^ P + ^3SSin2 + ^ X,){Ml^ - XpS^) 



+ 



2A3s(A„d+A^s) 
sin 2/3 



In 



I)} 



{m-^nf^ ^ 3^ { [^^%#^ (A.ssin^ /3 + Asscos^ + (A. + X,){Ml^ - X^s^) 



K) 

(Ml)" 



32fr2 I 



^Xp{Mly^ 



{A-uc 



Ar., 



X^X'pS sin2/3 
^167r2(A„d+Ajfs) 

167r2(A„rf+Ajfs) ^■^^ 

A4^ud /]„ / \ 
327r2 A4 / 



(Aw + 3A^s) In 

3A|f s) In 



sin 2;8 



sin 2,9 



t)} 

In 



A" 



ud 



In 



A4 ApS sin 2/3 
327r2(A„d + Ajfs) 
A,A^ssin2^* 



(.A^fj 



QX^s) In 



(^S 



Ac 



327r2(A„d + Ajfs) 
A3A^ssin2/3 
' 327r2(A„d+A^-fs) 
3A^ ll ■ 



(3Aw + 4A:^s) In 
(3Aw + 4A^s)ln 



327r2 \ A" 

9ApAyssin2ff 



In 



167r2(A„d + Ajfs) 



(A 



-F3A^s) In 



In 



In 



A* 



2A>s(A„d+A^s) , 
sin 2/3 
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1 

327r' 



^{(A2 + Ai)ln(^ 
'(A2 + A2)ln/'"'" 



(Ai)^ 

(A3 + A4)^^ ^^[{2XXs + Xl) 



(As)''^ ~ 3^ sin 2/32 
(Ae)''* ~ ||^sin2/31n 
(At)""* ~ ^ sin 2/3 In 




(A|>)^" 
(A|.)«" 



327r2 

32^|Ap(A.+A3)ln< 
Ap(A. + A3)ln( 



1 

327r2 



A" 



+ Ap(A. - A3)ln< 



+ Ap(A3 - A.)ln< 



{Xm)" ^- 



\\p sin 2/3 /Q 4 j+A^q^lIn 
X3X%' sin 2 f3 ,4-\M^\]„ 



A4 A}i, sin 2/3 
327r2(A„d+Ajfs) 



(2A„d + 3A^s) In 



327r2 



{ 



In 



A4 



In 



} 



Ap \ A4Ap A„dsin2;3 i 

(«2)H±^_L_{(AM2 + ^^2)i,^,« 

/'^2\^=^ A^Aj. A^dsin2;3 i„ / \ 

l%j - -8^^^(A„<j+A^«) ^^\^/ 

(4)''" ^ - lt-,^::f/.,, (4A„, + 3A^5)ln 



X^Hn 



Charged-Higgs contributions — nMSSM: 

The Charged-Higgs boson squared mass-matrix now reads : 

Hi,) = + (A^ds - ml,) tan"! /3 + AJF^p + X^' + Xpi\S\^ - s') 

{MlMS,Hl,) = Ml^ + {A^as - m?2) tan/3 + A.|iJ°|2 + Aal/f^p + Xp{\S\'' - s') 
{Mli)r2{S, H%) = Au,S - mf2 - MiHlH^af 

Applying the same recipe as in the previous paragraph, we obtain the corrections: 



327r2 
1 

3271-2 



- 32^^^{ 



H 



A" 



In 



'^^"'^'^""'^-^ (A. cos2 /3 + A3 sin^ /3) + (A. + A3)(M|^ - A^s^) 



sin 2/3 



2(A„da— 7»i2) 



2M{A^dS-ml2) 
sin 2/3 



In 



} 



(A. sin^ ^ + A3 cos2 p) + (A. + A3)(M|^ - ApS^) 



I 2A.(A„dS-mi2) 
sin 2/3 



A4m 



(ml) 

(mI) 

^ ^ {in 

/ . .... 



2Ai 



p 

mlm- 



A4 



In 



A4 



In 



Ap.s2 



I)} 



sin 2^ 



I)} 



In 



In 



A" 



A4 
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In 



A" 



, 2Ap(A„dS-7»i2) 1 

sin 2/3 



(C.25) 



Auc 

{A, 
(Ads) 
(As) 



[{As) ^ 



A4 Aj>AudS^ sin 2/3 
32Tr^ (^A^dS-mj^) 
^ A.A„dni?2 '5in2;g 

' ~ 327r2(A„dS-mf2) ^ 
ff± ^ A3^,.dm^2 sin 2/3 
~ 327r2(yl„,,«-m22) 
ff* _ 3Ap.-V,';rf,s- siii'^ 2^ , 
~ GATT^iA^dS-mj^)^ 

3ApA^dm^2 sin 2/3 
327r2(A„<js-"»f2) 
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(Ai) 

(A3 + Xi)" 



1 

327r2 



(A^ + Ai)ln 
32^UA2 + A2)ln 



(A, 
(Afi 

(At 



(A^)^" 
(A^)^" 



1 

327r2 

sin 2/3^ 



{(2A.A3 + Al)ln 




327r2 

Sin 2/3 In 

^sin2/31n(^ 
, 1 

■ 327r2 



Cxi. 



Ap(A. + A3)ln 
3^ |A>(A. + A3)ln 

(A3-A.)A„dAj,gsin^ 2/3 

647r2(/l„dS— "112) 
(A.-A3)A„dAj,ssin^2)g ^ 
647r2(A„dS-mj2) 

A4Apsin2/3 z^) < 
327r2(A„dS-m22) X'^^udS ' 
A4ApaA„d sin 2/3 



(Am) 

(\M\H^ 

yAp ) _ 327r2(A„dS-m^^ 




327r2 Y^P \ A4 
A^^A>ssin^2/3 , 
{AudS-ml^y 



mf2)ln 



Neutral-Higgs contributions - NMSSM: 

The 6x6 symmetric squared mass-matrix is given by its entries, in the base {h}l,h%h%,a^,a°^,a%): 

s 



{M 



2 ^ 



+ {A^d + XpS 



tan (3 



+ Xp{\S\' - s') + XX2\Hy + Rem) + (A3 + X,)\H'^ 



t0]2 



(A^ffo)22 
(^lfo)l2 

{M%o)55 
{Mjjo)45 
{M%o)66 
{M%o) 

(Mho) 56 

(A^ffo)24 
{Mho)25 

(Ml^o)34 
(A^^o)35 



Ml^ + (A„d + A|/5-)s tan/3 + Xp{\S\^ - s") + A.(2|i^;^|2 + Ro{H'}'^)) + (A3 + Xi)\Hl 



0|2 



-AwRe(-S) - A:^Re(52) + 2(A3 + A4)Re(ifO)Re(ff5) 

(^s + 4k2s)s + 2^s(Re(5) - s) + 2k'(2|S|2 + Re(5') - Ss^) + Ap(|if°|2 + - 2A^Re(F°ifS) 

-AwRe(ifS) - 2A^Re(S'*ilS) + 2A>Re(S')Re(iJ°) 
-A„dRe(if°) - 2A^Re(S'*i?°) + 2A>Re(S')Re(FS) 



M%^ + {Aud + X'^s) 



tan^ 



+ X-p{\S\' - s') + X.{2m' - R^iH^)) + (A3 + Xi)\H^^ 



7-0 12 



jfo;46 
2 



+ {A^a + A^s)stan/3 + Ap(|5|2 - s^) + A.(2|if0p - Re(/f02)) + (A3 + Xi)\H^S 
AudRe{S) + A^Re(52) + 2(A3 + X^)\m{H'i)\m{Hl) 

-^Ass - 2^s(Re(5) -8)+ 2k\2\S\^ - R^{S^) - s^) + Xp{\Hy + \Hl\'') + 2A^Re(if0ff0) 
AudRe{Hl) - 2X^Re{S*H^) + 2Xplm{S)lm{H^) 
AudRe{H°) - 2X^Re{S*H^) + 2XpIm{S)Im{H^) 

A.Im(ifr) (C.26) 

A„dIm(S') - A^^Im(S'2) + 2(A3 + Xi)Rc{H^)lm{H°) 
AudlHHd) + 2Af Im(5ffS*) + 2Xplm{S)Rc{H°) 
AudlMS) - X^'lmiS^) + 2(A3 + X4)Iin(H^)Re{H°) 
XMiHD 

A„dIm(F°) + 2X^lm{SH°*) + 2ApIm(5)Re(//°) 
Audim{H^) + 2X^lm{S*H^) + 2XpRe{S)lm{H°) 
Audlm{H°) + 2X^lm{S*H°) + 2A>Re(5)Im(ifS) 
-2AsIm{S) + 2K^lm{S^) - 2X^lm{H^H^) 
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One can expand its eigenvalues in terms of doublet fields: 

ml {S, Hi = m] (5) + {S, J?°) + {S, H°) + O {{H%f) 
The associated potential is then given by: 



(C.27) 



2(2), 



J-0 ttO\ 



A2 



In 



In 

3 

' 2 



A2 



2(1), 



J-0 TjO\ 



In 



A2 



A2 



The large logarithms arc then those terms multiplying In 



2 {0).„^ 
A2 



~ In 



2(0)/ \ 

~7P 



(C.28) 



+ . . . Thence consists 



our primary task in diagonalizing Mj^o perturbatively with respect to the doublet fields. 
Wc first consider M'^^^){S,H^^ = 0) in order to obtain mf^^\s). We denote as \Ei) = {Sij)j=i^,,,fi the 
elements of the canonical base of M®. The subspaces Vec{|£^i) , {£2} , I-E4) , I-E5)} and Vec{|£^3) , \Ee)} 
obviously decouple in M'j^oiS, H'^ ^ = 0). In the doublet sector, one notices that the eigenstate equation 

M]jo{S,H'^^^ = 0) - m^^ Z]i=i,2,4,5 = is equivalent to ^7W^ - m^^ (-^1,-22)^ = 0, where zi = 

xi + 1x3, Z2 = X2 — 1x4^ and AA^ is the 2x2 (complex) hermitian matrix determined by the following 
entries: 



tan /3 



+ \-p{\S\^ - s^) 



M'i2 = + {And + X^'s)s tan P + \p{\S{' - s 



(C.29) 



Mj2 



~-^udS 



One recognises A^^±(5, if^^ = 0), up to the sign of the off-diagonal terms. Ai'^ is diagonalized by the 
eigenstates: 



{zi)ho = 



m 



HO 



Ml, 



= Xd 



sin 2/3 



-j +(|Aw5*+A^52|'-(Ad + A^s)V) 



1/2 



-Ml,* 



,J{ml,-Ml,Y + \Ml^ 



2 12 



= -Vd 



Ml. 



= Vd ; {z2)h° = 



Ml, 



'{ml,,-Ml,f + \Ml,\^ 
The following relations will proove useful later: 

2 {Aud + \'^s)s 



\l{mlo-Ml,f + \Ml, 



Xd 



2 12 



Xd + \yD\ 



Xd - \yD\ 



tan 2/3 mla{S) - m]jo{S) 



_ AudS* + A|f 



(C.30) 



(C.31) 



Zi = {{zi)i^o, {z2)ho)^ and Z2 = ((zi)j/o, (^2)//")"^ eigenvectors of Ai"^ in the complex sense. In the 
real sense, iZi and 1Z2 form two other linearly-independant (and degenerate to Zi, Z2) eigenstates. We 
thus obtain the doublet eigenstates of M.'^jjo{S , ^ = 0): 



\h1) = XD \Ei) - ReiVD) I-E2) - luiiVD) I-E4) 
j/i^) = -lm(yo) |i;2) + XD 1^3) + Re(yB) IS4) 

= Re(yD) + xd 1^2) + Im(yzj) 1^3) 
\H^) = -lm(yc) + Re{yD) 1^3) - a;^ I-E4) 



(C.32) 
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For the remaining singlet-states, Ai^(,{S, ^ = 0) is diagonalized by: 

= i^As + 2k's)s + 4K'i\S\' _ s^) + / - 2 [\AsS + k'S*'\'] ^'^ 



\h%)=xs\E^)-ys\E^) 
V''c>i=ys\Ez)^xs\E^) 



(C.33) 



A^^o(5',0)66 - ■m\a 



ys 
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[M\o{S,QU-mloJ +\M\o{S,QU\^ ^ (x^„ [S, - K") ' + \^ho (S, 0)36 1 

It is convenient to introduce the notation zs = xs + lys and note the following relations: 



\zs\ 



1 



4 = 4- 



m{,{S)-mi,{S) 



(C.34) 



At this stage, one can already determine the pure-singlet pararneters of the potential. Similarly to the 



Rd 



1 



charged case (eq. C.21), one can formulate the first term in eq. 



m 



HO 



m 



hO 



J7 T5 = + m?o and Rs 



C.28 



in terms of To 



[m 



+ mlo), 



m 



m 



Moreover, the contributions from the 



doublet are trivially identical to those in the charged case and here, as well, we give only the leading 
term in sin 2/3 — )• for the coefficient multiplying In ( ). The logarithmicaly-enhanced parameters are 
then: 



(At) 



H" 



-3M:.::^(5A„, + m-.)in 

+ i [A| ~ lAsK^s ~ 4«:S2] In ( ^ 



Ass(As+1k''s) , 
12871-2 "1 



sin 2^ 



In 



2Aps(A„d+Ag's) 
sin 2/3 



In 



3A?. + 10A1.k^s+2A?,k*s2-60A2„k''s^-20Ask**s^+8k"'s^ 



3\p XpS-' sin2i3 



167r2 A„rf+A-J;'s 



i2Aud + SXjIs)\n 



3 As A%+2A% k'^ s-6A%K^ -32As ~16k^ 



In 



2567r2 



(4) 



s s 



A" 



16 



{As+K^s)3 



In 



In 



10247r2 ■ 



9ApApSsin2^ .^-^XAf^Nl^ 
16,r2(A„, + Aj.^.)^^"'^ + '^^-P 



9k^As A% + 17A%k^ s+4:2A%k* -8As k*^ -16k^ s'^ 



k'^s{As + k^s)^ 



In 



5 In 



rar.m 



^ 647r2 1 \ A" 

.H" ^ A^^Ap A„dsin2ff 



j_ A^s + liA%K^ s+?,2A%K.'^ +2b6A%K'^ -GiOAs K.'^ s'^ -2bGK^° s'^ , 

64 Kls2(yls+K;2s)3 



> M 



fAjj sin 2/3 



In 



"'"^'^ (4 A 



K^ylg 5At+22AiK''s + 144AiK*5- -352AsK s -128k'*s*' 



20487r2 



3A*^s) In 



In 



K-'As A|+2A|k2s-36A|k'*s2-344Ask''s^-64k**s* 



5127r2 



2(As+«;=s)3 



In 



(C.35) 



The next step consists in considering 0{H^) corrections to the neutral eigenvalues. Note that, for the 
eigenvalue m^^^^\s) of M.'^^q{S,H^^ = 0), with Tr, the trace on the corresponding eigenspace. 



(C.36) 



where M^j^l'"^^ stand for the matrices with terms of 0{H^''^) in TW^oj Pj corresponds to the projector on 
the eigenspace of the eigenvalue m'^^^\ Defining r] = 1; —1 for D = H^; and e = 1; — 1 for S = h^]a^ 
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and using the relations C.31 C.34 we obtain the following matrix elements: 
Tr,, {D\ Mlf \D) = (Ds) + ^ ^ {Da) 



'iff '"'ho 
{Ds) = [2X.+Xs + X,]i\H'j' + \HX) 



{S\M'^^^\S)^{Ss) + 



2{Aud + \Vs)s 
tan 2/3 



[2X.-X^-Xi]{\Hl\ 
(Sa) 



- 4(A3 + A4) [A^dRe{SHliH°) + Xl!Re{S*^H°H° 



a% 



(5s)=Ap(|i7°p + |i7°n 

{Ss) - -8A^;^ [AsRe{SHlH',) + ^^Re{S* ' H^H^,)] 
{S\m]}o''PdM 



2(1) D„ W \S) ^ (1; 1) + ^ (-1; 1) + — 

■ 2| c|21 n rr0|2 I I rr0|2\ 



(-i;-i> 



(1; 1) = 2 [Ala + 4A^^2|5p + 2Ap2|5p} {\Hl\^ + - 2Ap [A„,Re(5i7„"ff°) + 2A^^Re(5* ^iJ^ff")] 

(-1; 1) = 2Ap [A^l^p + 2X'JM^ + iX^^A^dRe{S^)] {\Hy + \H^,\^) (C.37) 

- 4Ap2|5p [A„,Re(5i7°i/S) + X'jRe{S*^HlH°)\ - 8A|/A,d [A*^|5|2Re(5i/°iJ,") + A.^dRe{S* ^ H^H^)] 
(1; -1) = [8A*f A„d(As|5|2 + K'ReiS^)) + 4Ap2 (AsRe(53) + K'\St)] {\Hy + 

- 8Ap [{K^Aud + 2A|f As)|5|2Re(^i/077S) + {AgA^a + 2A^^K2|^|2)j^g(^* 2^o^O)] 
(-1; -1) = 8Ap [AU^^'^A^d + i\¥As)\S\^ + 2X^J^t^\S\^ 

HAsAl, + 3A|,^A„,Ac2|5p + 2A|f 2^s|5nRe(53)] (IH^P + W'd?) 



S{Aud + A|f g)g 
tan 2/3 



16Ap2 (AsRe(^3) + K^l^n [A„dRe(5ffOi/0) + A|fRe(5* ^i/OiJ^)] 



8 [{X'MsiAld + AX^'W) + n'Al,) |5pRe(5i/°i/°^ 

+A^d {AsiAl, + 4Xj^'\S\') + AXj!^^'\S\') Re(5*2i7,«i/°) 

+Xl!Al,n'Re{S'H°H',)+4Xl!'A^dn'ReiS*'H:H°) 



Getting back to C.28 one obtains the 0{H'^) coefficients: 



64^ { [Aid + 4A. (m|>^ + Aps2 + 2 cos^ (3 



2 fl (A„d + Ag's)s 
sin 2/3 



+2(A3 + A4) (m|^ - Aps2 + 2 sin2 (^(A^^+^l^ 



sin 2^ 



A* 



+ [Al^ - 2Aps(A5 + 2k^s)] In 
6^ { [Aid + 4A. (m|^ + Aps2 + 2 sin^ /3 



2 o(A„d + Aj.-'s)s 



+2(A3 + A4) (m 
[AL - 2Aps(As + 2k^s)] In 



2 



sin 20 

„2 I 0„^o2 fl(^..d + Ag.-'a)s 



Aps^ + 2 008-^/3^ 



sin2^( 



A* 



In 



A* 



In 



A* 



{A^^df ~ 32^|(A3 + A4 + Ap)A„<iln(^^) + (ApA„rf + 2AffAs)ln/^^ 
(A^)^'' ^ |(Ap(2A.A3 + A4 + Ap) + X'^P) In (^) + 2 (Ap(4.j2 + Ap) + 2A|.^2) 
(Al,)^'' ^ |(Ap(2A.A3 + A4 + Ap) + XP) In ) + 2 (Ap(4.j2 + + ^yM2^ 

iXp"" ^ |a|.^(A3 + A4 + 2Ap) In ) + 2Xj.\K^ + Ap) In ^^^^i^ + 



m "I n 

In/ "s °g 



In 



A4 
2 2 

A^ 



(C.38) 
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The logarithms involving ratios of Higgs masses (symbolized by . . . ) are too complicated to write 
down explicitely. They also appear within contributions to the Zs-violating parameters. We could 
check however that such contributions to Zs-violating parameters vanished in relevant limits {rn^o ^ 

m|o sin2;3 0, m^o < ^h%,a%' '^1% ~^ 

We skip the computation of corrections to the quartic doublet parameters, here as well as for the neu- 
tralino contributions: such a task, although straightforward (perturbative calculation of the eigenvalues 
of a matrix up to the fourth order) promises to be technically tedious. 

Neutral-Higgs contributions - nMSSM: 
Wc can draw some conclusions from the study in the NMSSM case. For simplicity, we will confine here 
to contributions to pure singlet parameters, so that we need only consider Ai^jjo{S, H^^ = 0). Obviously, 
doublet and singlet sectors will decouple again, and the doublet sector will generate the same corrections 
to the Aj's as the charged Higgs sector. As for the singlet sector, it is particularly simple in the nMSSM 
since no dependence in S appears. We thus obtain the leading-logarithms to the pure-singlet coefficients: 




(C.39) 



SumniEiry of the analysis: 

2 

Among the potentially large logarithms, we may distinguish among those of the form In^, which 
compare a given sector to the scale A, typically chosen as the mass of the third-generation squarks, or to 
another sector, and those sensitive to hierarchies within a sector In — I-. 

2 

In the case of the NMSSM, the logarithms In ^ obviously appear only in the corrections to the 
Za-conserving parameters. Moreover, when logarithms of the type In — ^ appear within Za-violating 

9 9 
171 ■ TTl- 

parameters, they tend to be balanced by prefactors vanishing in the hierarchical limit (typically (^2'_^2^ ); 

so that they cannot be regarded as an enhancement factor (contrarily to when they appear in Z3- 
conserving parameters, where the prefactor does not necessarily vanish in this limit). One can thus 
conclude that leading-logarithms preserve the Zs-induced structure of the potential. 

For the nMSSM, the Za-symmetry is actually still present at tree-level in all the sectors of the 
spectrum, with the exception of the Higgs sector, where it is explicitly violated. Consequently, large 
logarithms still favour the Za-conserving terms (even though they are not all present at the classical 
level), while Za-violating effects perdure in the Higgs sector. In that case, large logarithms seem likely 
to destroy the classical structure. 
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